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Primary  activities  during  the  second  year  consisted  of  (i) 
expanding  the  laboratory  capabilities  for  processing  and  testing  of 
composites,  (ii)  conducting  research  in  accordance  with  the  Statement 
of  Work  given  in  Section  1.3,  (iii)  preparing  nine  technical  papers 
and  M.S.  theses,  and  (iv)  various  interactions  of  the  faculty  with  the 
technical  community  through  presentation  of  papers,  participation  as 
members  of  technical  committees,  etc. 

Sections  2-5  summarize  the  research  activities.  The  professional 
personnel  associated  with  the  project  and  the  outside  activities  of  the 
faculty  related  to  composites  are  given  in  Section  6.  Papers  completed 
and/or  published  during  the  year  are  reproduced  in  Appendix  A.  A 
brochure  describing  the  student  activities  is  in  Appendix  B. 

1.2  Discussion 

The  general  objective  of  the  research  program  is  to  develop  improved 
understanding  of  structural  composite  material  as  a  basis  for  improving 
its  structural  performance.  Integral  to  achievement  of  this  objective 
is  the  determination  of  constitutive  equations  and  time  (or  cycles)  to 
fracture  for  complex  loading  patterns  and  hostile  environments  such  as 
high  relative  humidity.  This  in  turn  will  require  a  better  understanding 
of  the  viscoelastic  deformation  and  micromechanisms  of  fracture  in 
composites  including  delamination,  cumulative  damage  via  microcrack 
Initiation  and  growth,  and  meaningful  failure  criteria.  The  importance 
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of  micromechanisms  of  viscoelastic  deformation  and  fracture  in  epoxy 
materials  has  only  recently  become  widely  appreciated;  for  example, 
graphite/epoxy  composites  are  quite  sensitive  to  processing  and  to 
service  environments,  especially  moisture  [1].  Furthermore,  because 
curing  is  seldom  taken  to  completion  in  the  epoxies  used  in  graphite/ 
epoxy  composites  [2],  the  degree  of  cure  is  a  processing  variable. 

Physical  aging  of  the  epoxy  matrix  may  occur  in  service,  giving 
significant  changes  in  mechanical  properties  [3].  The  post-cure 
cool-down  path  affects  the  magnitude  of  residual  stresses  in  a  composite 
[4].  Moisture  absorption  can  change  the  glass  transition  temperature, 
giving  a  significant  change  in  mechanical  properties  [1].  The  rate  of 
moisture  absorption  appears  to  depend  on  the  degree  of  curing  and 
physical  aging  (densification) ,  and  moisture  content  may  itself  affect 
the  rate  of  physical  aging. 

The  lack  of  consistency  in  experimental  results  and  micromechanistic 
modelling  in  the  epoxy  matrix  composite  materials  literature  today  is  a 
consequence  of  the  fact  that  much  of  the  earlier  research  was  conducted 
without  fully  controlling  some  of  the  more  subtle  variables.  It  is 
anticipated  that  the  continuing  research  program  at  Texas  A&M  University 
described  herein  will  make  a  major  contribution  toward  a  more  unified 
understanding  of  composite  materials  behavior. 

In  view  of  the  large  number  of  variables  to  be  considered,  it  is 
essential  that  the  mechanical  testing  be  guided  by  mechanistic 
modelling  to  keep  the  total  amount  of  experimental  work  at  a  manageable 
level.  A  considerable  part  of  the  effort  at  Texas  A&M  University  is  being 
so  directed. 
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1 . 3  Statement  of  Work 


•  "a.  Investigate  the  effects  of  composite  cure  cycle  parameters  on  the 
mechanical  characteristics  of  structural  specimens  and  components  (see 
Section  2. ): 

(1)  Investigate  the  development  of  new  innovative  in-process 
monitoring  devices  and  study  improvements  in  the  "Duomorph"  sensor. 

(2)  Conduct  experimental  measurements  on  both  neat  resin, 
unidirectional  laminates  and  cross-ply  laminates  to  establish  the 
relationships  at  various  elevated  temperature  and  relative  humidities 
between  processing  parameters  and  mechanical  strain  rate  behavior, 
linear  and  identified  nonlinear  viscoelastic  behavior,  dimensional 
stability,  and  fracture  characterization. 

(3)  Investigate  the  influence  of  processing  parameters  on  T 

at  various  relative  humidities.  9 

•  b.  Investigate  the  influence  of  elevated  temperature  and  high 
relative  humidity  environments  on  composite  laminates  and  adhesively 
bonded  structural  components  (see  Section  3.): 

(1)  Investigate,  on  a  theoretical  level,  the  formulations  of  a 
coupled  diffusion  theory  to  accurately  predict  the  environmental 
degradation  of  composite  laminates  and  adhesive  mechanical  performance 
which  considers  the  influence  of  cure  cycle  parameters  being 
investigated  in  a  above. 

(2)  Predict  moisture  concentration  profiles  and  internal  stress/ 
strain  fields  within  a  laminate  and  bonded  structure  as  well  as  in 
the  vicinity  of  exposed  edges.  These  profiles  will  consider  generic 
environmental  conditions  of  absorption  and  desorption. 

(3)  Compare  theoretical  predictions  with  available  experimental 
data  or  define  experiments  which  would  yield  necessary  data  for  this 
comparison. 

•  c.  Investigate  suitable  experimental  techniques  which  can  be  used  to 
nondestructively  examine  composite  laminate  initial  defects,  damage 
growth  and  environmental  degradation  characteristics  which  affect 
mechanical  performance.  Initially,  emphasis  will  be  given  to  neutron 
radiography  as  a  suitable  technique  (see  Section  4): 

(1)  Conduct  necessary  experiments  to  characterize  neutron  beam 
parameters  and  associated  parameters'  sensitivities  applicable  for 
nondestructively  examining  major  composite  laminate  mechanical 
characteristics  which  affect  strength,  dynamic  response  and  durability. 

(2)  Investigate  the  development  of  experimental  neutron 
radiography  techniques  to  measure  moisture  concentrations  in  composite 
laminates  as  well  as  spatial  concentration  profiles. 


(3)  Investigate  the  development  of  automated  procedures  for 
the  analysis  of  data  from  neutron  radiography  experiments  which 
accurately  and  efficiently  characterize  and  identify  laminate 
initial  defects,  mechanical  damage,  and  moisture  concentrations. 

(4)  Investigation  of  automated  data  analysis  procedures  in  c(3) 
above  will  consider  both  planar  and  stereographic  data  displays." 


jk 
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2.  PROCESSING  AND  TESTING  OF  RESIN  AND  COMPOSITES 

2.1  New  Equipment 

The  processing  and  testing  capability  is  being  expanded  through 
the  acquisition  of  a  large  amount  of  new  equipment.  Recently  acquired 
is  a  MINC-11  laboratory  data  acquisition  and  control  microcomputer 
system  (manufactured  by  Digital  Equipment  Corporation)  and  a  digital 
plotter  to  provide  high  quality,  permanent  graphic  output.  Orders 
have  been  placed  for  1)  curing  oven,  2)  vacuum  oven,  3)  automatic  controls 
for  use  on  existing  press/autoclave,  4)  cold  stage  for  differential 
scanning  calorimeter,  5)  diamond  tools  for  machining  composites, 

6)  humidity  monitor,  7)  digital  thermometer,  and  8)  laboratory  oscilli- 
scope.  Funds  were  provided  by  AFOSR  and  Texas  A&M  University. 

A  considerable  amount  of  additional  testing  and  data  acquisition 
equipment  is  to  be  purchased  during  the  coming  year  including  a  closed- 
loop  MTS  mechanical  testing  system,  two  data  acquisition  and  micro¬ 
computer  systems,  digital  oscilliscope  and  an  ultrasonic  tester  with 
C-scan  capability.  Funds  are  to  be  provided  by  Texas  A&M  University. 

2.2  Processing  of  Composites 

The  laminating  press/autoclave  fabricated  under  the  project  and 
placed  into  operation  last  year  has  been  used  extensively.  Faculty  and 

i 

student  research  assistants  have  learned  the  detailed  procedure  for  i 

1 

laminating  and  curing  a  structural  composite.  They  have  used  the 
press  to  fabricate  specimens  used  in  the  several  research  tasks. 


The  laboratory  data  acquisition  system  has  been  used  to  gather  data 
from  creep  and  recovery  experiments.  Data  analysis  software  employing  a 
nonl inear  Gauss  iteration  technique  for  estimating  parameters  in  the 
following  representations  for  creep  compliance  was  developed: 

0(t)  =  D 0  +  D-j tn,  D(t)  =  Dq  +  D^log  t,  and  D(t)  =  D^tn 

Many  plotting  routines  have  been  developed  which  produce  high  quality 
and  high  resolution  graphic  presentations  of  data.  Figure  1  is  an  example. 


2.4  Dynamic  Mechanical  Response  of  Resins  and  Composites 

This  portion  of  the  study  on  the  effect  of  cure  cycle  parameters  is 
concerned  with  the  direct  monitoring  of  the  viscoelastic  complex  modulus 
of  resins  and  composites  during  and  immediately  following  the  cure  cycle. 

We  have  adapted  a  piezoelectric  device  called  the  "duomorph"  for  such 
use. 

The  duomorph  consists  of  two  crystalline  piezoelectric  plates  which 
are  bonded  together;  a  central  metallic  plate  is  sometimes  used  to  alter 
the  response  characteristics  of  the  assembly.  The  duomorph  is  electrically 
excited  and,  through  a  study  of  its  bending  response  when  in  contact  with 
a  material  of  interest,  the  viscoelastic  properties  of  this  material  are 
obtained.  The  gage  was  originally  developed  for  the  purpose  of  measuring 
in  situ  complex  modulus  when  embedded  in  the  fuel  in  a  solid  propellant 
rocket  [5];  for  this  application,  which  involved  telemetry  to  trigger  and 
measure  the  gage  response,  one  piezoelectric  plate  was  employed  as  a 


Figure  1:  Creep  strain  as  a  function  of  time. 


driver  and  the  voltage  response  of  the  other  was  used  to  obtain  the 
strain  response. 

The  latest  version  of  the  duomorph  developed  on  the  project  has  been 
demonstrated  to  perform  reproducibly  in  air  over  the  temperature  range  of 
25  tol80°C,the  upper  limit  being  the  approximate  cure  temperature  for 
advanced  composites. 

It  consists  of  two  1  inch  diameter  x  0.010  thick  piezoelectric  discs 
bonded  together  using  a  conductive  epoxy.  Evaluation  of  piezoelectric 
material  from  different  suppliers  has  resulted  in  the  selection  of  a 
new  material,  K-350,  which  is  a  lead  zirconate  titanate  piezoelectric 
ceramic  manufactured  by  Keramos,  Inc.,  Lizton,  Indiana.  As  with  the 
gage  used  in  the  study  of  propellant,  one  plate  serves  as  a  driver  and 
one  as  a  response  transducer.  However,  in  order  to  obtain  complex 
moduli  with  sufficient  accuracy  for  use  with  stress  analysis,  we  have 
found  it  necessary  to  work  with  the  ’iput  and  response  current,  rather 
than  voltage,  in  deriving  complex  moduli. 

Tests  of  Hercules  3501-6  epoxy  and  a  room- temperature  cure  epoxy 
resin  were  sufficiently  successful  to  justify  further  refinements  by 
upgrading  the  input/response  electronics  and  the  addition  of  an  automated 
data  acquisition  system  (which  will  include  a  Fast  fourier  transform 
analysis  capability).  Necessary  equipment  will  be  ordered  during  the 
next  quarter.  Also,  it  is  planned  to  adapt  the  gage  for  use  with 
laminates  and  develop  the  necessary  additional  theoretical  analysis  and 
experimental  techniques.  The  high  stiffness  of  composites  compared  to 
the  neat  epoxy  will  probably  require  some  new  piezoelectric  (PZT) 
elements  and  gage  design.  One  possibility  would  be  a  PZT  gage 


mounted  in  a  recessed  cavity  in  the  bottom  metal  plate  which  supports 
the  laminate  during  cure.  The  gage  would  rest  on  a  relatively  soft  layer 
so  that  the  only  resistance  to  stretching  and  bending  would  arise  from 
contact  with  the  laminate.  The  PZT  material  would  be  polarized  so  as 
to  deform  primarily  in  one  direction,  and  the  gage  would  be  placed  so 
that  it  is  at  90°  to  the  fiber  direction  of  the  bottom  ply.  Through 
analysis  and  calibration,  the  amount  of  bending  and  stretching  will  be 
used  to  extract  the  insitu  complex  modulus  of  the  epoxy,  or  at  least  a 
relative  measure  of  this  modulus.  It  is  expected  that  before  and  during 
gelation,-  the  changes  in  bending  resistance  will  be  most  sensitive  to 
resin  changes.  However,  little  bending  will  occur  after  the  resin  is 
very  stiff,  and  then  the  amount  of  stretching  will  be  the  best  indicator 
of  resin  change.  Initially,  unidirectional  layups  will  be  employed. 

After  the  method  is  developed  for  this  special  case,  laminates  with 
various  ply  orientations  will  be  monitored;  with  the  intent  of  using  the 
gage  response  in  a  closed-loop  control  of  the  press/autoclave. 

2.5  Physical  Aging  of  Epoxy  Resins 


The  study  of  aging  phenomena  following  curing  and  quenching  from  the 
glass-transition  temperature  T^  is  continuing.  The  resin  matrix  of  the 
composite  is  below  the  glass  transition  in  use  and  so  continues  to  change 
properties  with  time  due  to  a  relaxation  of  the  free  volume  toward  its 
equilibrium  value.  These  changes  can  be  observed  in  a  variety  of  ways; 
we  are  doing  It  primarily  by  measuring  the  specific  heat  curve  near 
Tg  in  the  differential  scanning  calorimeter  (DSC)  and  by  studying  the 
creep  and  recovery  behavior  of  flexure  samples  in  the  thermo-mechanical 
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analyzer  (TMA).  Some  preliminary  investigations  with  the  duomorph 
(cf.  Section  2.4)  indicate  that  this  device  will  prove  useful  in  the 
study  of  aging,  and  therefore  this  approach  will  be  further  explored 
during  the  coming  year. 

The  changes  in  the  specific  heat  which  we  observe  are  quantities 
which  can  be  related  to  several  molecular  models  of  the  glass  transition, 
which  imply  that  the  mechanical  characteristics  such  as  brittleness  will 
also  change  in  time.  Such  changes  would  be  of  importance  in  composite 
materials  since  aging  has  been  found  to  continue  over  many  decades  of 
time  in  other  thermosetting  materials.  Any  mechanical  characteristics 
of  the  composite  which  depend  upon  the  brittleness  of  the  resin  can  be 
expected  to  evolve  over  a  matter  of  years  if  our  suppositions  are  correct. 
A  variety  of  resins  have  been  studied,  which  includes  the  proprietary 
resins  as  well  as  a  number  where  the  crosslink  density  can  be  estimated. 

The  work  is  tedious  since  moisture  content,  completeness  of  cure,  and 
other  factors  also  influence  the  same  experimental  quantities.  However, 
the  work  is  sufficiently  far  along  that  we  expect  to  be  able  to  prepare 
a  technical  report  on  the  DSC  studies  during  the  next  quarter. 

The  aging  studies  conducted  with  the  TMA  are  described  in  the  M.S. 
Thesis  "Determination  of  the  Relationship  of  Free  Volume  to  Mechanical 
Behavior  for  an  Epoxy  System  Subjected  to  Various  Aging  Histories" 

(cf.  Section  5.2). 

2.6  Mechanical  and  Failure  Properties  of  Composites  Including 

gnv^ronmentai  Effects 

Various  investigations  have  been  completed  and  are  reported  in  M.S. 
Theses  (cf.  Section  5.2)  and  in  the  paper  "A  Method  for  Determining  the 


Mode  I  Delamination  Fracture  Toughness  of  Elastic  and  Viscoelastic 
Composite  Materials",  which  appears  in  Appendix  A. 

Further  work  on  mechanical  and  failure  properties  will  be  conducted 
during  the  coming  year,  as  indicated  in  the  list  of  M.S.  Thesis  topics  in 
Section  5.1. 


3.  ANALYSIS  OF  COMPOSITES 


3.1  Summary 

Many  analytical  results  have  been  obtained  for  elastic  and 
viscoelastic  stresses  and  deformations  due  to  moisture  and  temperature 
changes  in  composite  materials  and  for  optimum  cooling  paths  following 
cure;  they  are  covered  in  papers  given  in  Appendix  A  and  the  M.S.  Thesis 
"Stresses  due  to  Environmental  Conditioning  of  Cross-Ply  Graphite/Epoxy 
Laminates"  (cf.  Section  5.2).  Energy  release  rate  analysis  of  nonlinear 
cantilever  beams  was  conducted  in  support  of  the  delamination  fracture 
toughness  studies  (cf.  Section  2.6). 

During  the  coming  year  elastic  and  viscoelastic  stress  analyses  and 
fracture  analysis  of  composites  will  be  conducted  in  combined  theoretical  - 
experimental  investigations.  The  experimental  work  will  serve  to  confirm 
theoretical  predictions  and  guide  the  development  of  theoretical  models 
for  nonlinear  viscoelastic  mechanical  behavior,  physical  aging, 
environmental  effects,  and  damage  growth  and  fracture.  Student  M.S. 

Thesis  topics  listed  in  Section  5.1  indicate  the  primary  investigations 
to  be  undertaken. 

In  addition  to  analytical  work,  a  finite  element  model  for  the 
prediction  of  nonlinear  viscoelastic  response  of  polymeric  composites 
and  other  materials  is  under  development,  as  discussed  below. 

3.2  Nonlinear  Finite  Element  Model 

A  nonlinear  finite  element  model  [6]  was  generalized  in  the  previous 
reporting  period  to  nonlinear  viscoelastic  behavior  using  the  two- 
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dimensional  plane  stress/plane  strain/axi symmetric  case  and  programmed 
into  the  AGGIE  I  computer  program  [7].  The  current  version  of  the 
model  allows  several  alternate  input  options.  The  material  model  may 
be  either  nonlinear  viscoelastic  or  linear  viscoelastic  and  includes  as  a 
special  case  linear  thermoreologically  simple  behavior.  The  input  for 

r 

the  nonlinear  material  parameters  g^(o),  g2(a)  and  ao(a)  (see 
[6]  for  their  description)  may  be  in  terms  of  either  tabular  data  or 
equation  form  as  a  function  of  stress.  The  computer  coding  has  been 
checked  thus  far  by  running  two-dimensional  cases  which  degenerate  to 
a  uniaxial  stress  state. 

Further  work  will  include  verification  of  the  model  for  two-dimensional 
stress  states  and  extension  of  the  model  so  that  material  properties  are 
functions  of  other  parameters,  e.g.,  temperature  and  moisture.  The  AGGIE 
program  allows  for  general  continuum  analysis  through  the  use  of  4  to  8 
isoparametric  elements.  Hence,  the  program  will  allow  the  solution  of 
many  geometries,  including  laminated  composites.  Future  work  will 
consider  such  applications.  It  is  also  intended  that  the  nonlinear 
viscoelastic  properties  obtained  in  experimental  investigations  during 
the  coming  year  will  be  used  with  this  finite  element  code  to  predict 
response  of  composites  to  mechanical  and  environmental  inputs. 
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4.  NEUTRON  RADIOGRAPHY 


4.1  Summary 

Last  year  the  Nuclear  Science  Center  developed  an  extremely  high 
quality  neutron  radiography  facility  at  the  research  nuclear  reactor 
at  Texas  A&M  University.  The  operational  parameters  such  as  geometric 
resolution,  beam  uniformity,  film  density,  and  neutron  flux  were 
optimized  to  a  point  where  variation  in  material  density  of  less  than 
0.5%  may  be  measured.  A  variety  of  laminate  sections  with  known 
moisture  content  were  radiographed  in  an  attempt  to  establish  a 
relationship  between  film  density  changes  for  increased  moisture  in 
the  samples.  A  laminate  with  known  defects  was  also  studied. 

This  year  studies  were  directed  toward  the  improvement  and 
standardization  of  the  sensitivity  of  neutron  radiography  to  absorbed 
water  in  laminates.  By  using  polyethylene  to  simulate  water,  it  was 
estimated  that  at  least  1%  absorbed  water  was  required  for  detection. 
Subsequent  studies  with  laminates  saturated  to  this  level  indicated 
that  an  even  greater  amount  of  moisture  was  needed.  Inasmuch  as  the 
moisture  level  requirements  for  detection  of  moisture  gradients  are 
well  above  what  can  be  expected  in  laminates  in  realistic  environments, 
it  was  decided  to  terminate  this  area  of  investigation.  However, 
because  metals  are  transparent  to  the  neutron  beam  while  materials  with 
hydrogen  are  not,  this  facility  should  be  very  useful  for  flaw 
detection  in  resin  and  fiber  reinforced  resins,  especially  when  other 
detection  methods  are  made  ineffective  as  a  result  of  metal  barriers. 
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4.2  Beam  and  Material  Characteristics 

The  final  values  of  beam  parameters  that  have  been  reached  as  a 
result  of  the  research  effort  are: 

7  2 

Thermal  neutron  flux  >10  n/cm  /sec 

Scattered  and  fast  neutron 
content  1 0% 

6  2 

Gamma  intensity  4.0x10  n/cm 

mr 

Beam  uniformity  +3.5% 

Beam  divergence  1.5° 

Also,  optimum  film  developing,  handling  and  inspection  procedures 
necessary  for  reproducible  data  collections  have  been  established. 

Effective  neutron  cross  sections  for  all  of  the  materials  studied 
were  measured  and  are  listed  below: 


Material 

Et(cm"^ ) 

Scotch  ply  (GL/E) 

1.178 

Neat  Resin  (1) 

1.106 

(2) 

0.648 

Defect  Test  Plate  (GR/E) 

0.179 

Polyethylene 

1.084 

Water 

0.745 

(It  should  be  noted  that  these  cross  sections  are  valid  for  the  present 
geometry  only,  l.e. ,  relatively  thin  plate  samples  in  contact  with  the 
film  cassette.)  Also,  the  relationship  between  film  density  and 
integrated  neutron  flux  has  been  established: 
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D  =  e  (-0.6035  *R  +  0.1561  *R  -  0.0758) 
where  D  =  film  density 
$R  =  Relative  fluence 

Note  that  for  a  sample  of  thickness  t  exposed  to  a  neutron  fluence 

o  t 

*R,  the  transmitted  fluence  4>R  may  be  determined  from 


Thus,  the  neutron  cross  sections  and  these  two  equations  may  be  used  to 
determine  changes  in  sample  thickness  (At)  from  changes  in  film  densities. 
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5.  GRADUATE  RESEARCH  ASSISTANT  ACTIVITIES 


The  first  group  of  graduate  engineering  students  to  participate 
in  the  current  research  project  entered  the  program  in  September,  1978 
and  completed  the  requirements  for  a  degree  of  Master  of  Science  in 
August  1979.  Results  of  their  research  are  reported  in  the  following 
theses : 

1.  Devitt,  F.  D. ,  "Delamination  Fracture  Toughness  of  a  Unidirectional 
Composite. " 

2.  Douglass,  D.  A.,  "Stresses  due  to  Environmental  Conditioning  of 
Crossply  Graphite/Epoxy  Laminates." 

3.  Grieger,  K.  A.,  "Experimental  Investigation  of  Moisture  and 
Temperature  Conditioning  of  C600/5208  Graphite/Epoxy  Composite 
Material." 

4.  Ring,  D.  S.,  "Determination  of  the  Relationship  of  Free  Volume 
to  Mechanical  Behavior  for  an  Epoxy  System  Subjected  to  Various 
Aging  Histories." 

5.  Rogers,  J.  B.,  "Aeroelastic  Tailoring  of  Composite  Materials." 


Abstracts  are  given  in  Section  5.2.  Copies  of  the  theses  will  be 
provided  upon  written  request  to  the  Principal  Investigator  (R.  A.  Schapery). 

The  second  group  consists  of  eight  (8)  students  that  started  work  on 
an  M.S.  degree  in  engineering  in  September  1979.  The  thesis  topics  are 
listed  below;  most  studies  involve  both  experimental  and  theoretical  work. 

1.  Investigation  of  residual  deformations  and  curvatures  in 

unbalanced  laminates  due  to  cool-down  from  cure.  Viscoelastic 
effects  are  included. 
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2.  Investigation  of  residual  deformations  and  curvatures  in 
unbalanced  laminates  due  to  moisture  sorption.  Viscoelastic 
effects  are  included. 

3.  Nonlinear  viscoelastic  characterization  of  a  graphite/ epoxy 
composite. 

4.  Characterization  of  delamination  fracture  toughness  and  its 
relation  to  microstructure. 

5.  Delamination  fracture  due  to  compressive  loading. 

6.  Development  of  microcracks  in  cross-ply  laminates  due  to 
mechanical  loading. 

7.  Prediction  of  microcracking  and  delamination  in  composites 
by  the  finite  element  method,  including  a  study  of  crack 
stability. 

8.  Physical  aging  of  epoxy  resin  following  cooling  from  the  glass 
transition  temperature. 

The  next  quarterly  report  will  contain  thesis  proposals  describing 
the  research  planned  under  each  of  these  topics.  Besides  conducting  this 
research,  the  students  are  involved  in  academic  courses,  as  described 
in  the  brochure  in  Appendix  B. 

In  addition  to  graduate  students,  a  few  select  undergraduate 
students  assist  in  the  various  research  activities.  The  participation  of 
these  students  aids  in  the  research  and  helps  to  acquaint  them  with 
composite  materials  prior  to  enrolling  in  the  graduate  program. 

5.2  Abstracts  of  M.S,  Theses  Completed  in  August  1979 
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ABSTRACT 

Oelamination  Fracture  Toughness  of  a  Unidirectional 
Composite.  (August  1979) 

Daniel  Frank  Devitt,  B.S.  A.&A.E.,  Purdue  University; 

Chairman  of  Advisory  Committee:  Dr.  W.  L.  Bradley 

An  analytical  technique  for  reducing  experimental  data  is 
developed  for  a  unidirectional  composite  in  order  to  determine  the 
delamination  fracture  toughness  for  the  crack  opening  mode.  A 
unique  characteristic  of  this  data  reduction  technique  is  that  the 
crack  opening  mode  is  isolated  from  the  other  fracture  processes. 
Isolation  of  this  fracture  mode  allows  for  the  direct  determination 
of  the  critical  energy  release  rate  (Gc)  and  the  critical  stress 
intensity  factor  (KIC). 

The  values  for  the  critical  energy  release  rate  and  the  critical 
stress  intensity  factor  were  calculated  using  linear  elastic  fracture 
mechanics  coupled  with  nonlinear  beam  theory.  The  critical  energy 
release  rate  is  found  from  utilizing  linear  elastic  fracture  mechanics 
as  it  is  related  to  the  release  of  the  stored  strain  energy  in  the 
system.  Nonlinear  beam  theory  is  required  to  calculate  the  strain 
energy  in  the  system,  since  large  deflections  and  rotations  are 
produced  experimentally  and  are  not  accounted  for  in  linear  beam 
theory. 
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ABSTRACT 

Stresses  Due  to  Environmental  Conditioning 
of  Cross-ply  Graphite/Epoxy  Laminates.  (August  1979) 

David  Alan  Douglass,  B.S.  A.&A.E.,  Purdue  University: 

Chairman  of  Advisory  Committee:  Dr.  Y.  Weitsman 

This  thesis  concerns  the  internal  residual  stresses  due  to  the 
environmental  conditioning  of  cross-ply  graphite/epoxy  laminates,  where 
moisture  is  induced  into  the  material  by  exposure  to  high  relative 
humidity  at  elevated  temperatures. 

The  stress  fields  resulting  from  conditioning  at  two  temperature 
levels  are  determined  and  compared  for  various  symmetric  cross-ply 
graphite/epoxy  laminates  by  means  of  both  elastic  and  viscoelastic  analy¬ 
ses.  The  formulation  considers  temperature-dependent  moisture  diffusion 
and  time,  temperature  and  moisture  dependent  stress  relaxation.  The 
computations  are  based  upon  recent  data  and  employ  realistic  values  of 
material  parameters. 

The  analysis  shows  that  the  viscoelastic  stresses  are  much  smaller 
than  those  predicted  by  elastic  analysis,  and  that  conditioning  at 
150°F  results  in  residual  stresses  which  are  smaller  than  those  reached 
during  conditioning  at  180°F. 


''  mm.  y? Vl  t 
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ABSTRACT 

Experimental  Investigation  of  Moisture  and  Temperature  Conditioning  of 
C600/5208  Graphite/Epoxy  Composite  Material'  (December  1979) 
Kenneth  Allen  Grieger,  B.S.A.E.,  Purdue  Univeristy; 

Chairman  of  Advisory  Committee:  Dr.  R.  A.  Schapery 

There  are  three  main  questions  to  be  investigated  in  this  thesis. 
This  first  objective  is  to  determine  if  a  very  high  relative  humidity 
(e.g.  98  percent)  environment  has  approximately  the  same  effect  as  the 
more  easily  conducted  water  soak  condition  on  the  moisture  gain  and 
damage  level  of  the  specimen.  The  second  objective  is  to  determine  if 
the  maximum  moisture  content  at  different  relative  humidities  is  inde¬ 
pendent  of  tempeature.  The  third  objective  is  to  determine  the  amount 
of  surface  cracking  for  various  environments. 

An  examination  by  means  of  a  light  microscope  and  water  absorption 
rate  of  C600/5208  composite  laminates,  both  crossply  and  unidirec¬ 
tional,  which  had  been  exposed  to  high  temperature  and  moisture  con¬ 
ditions  is  described.  An  attempt  is  made  to  explain  these  observations 
by  means  of  surface  cracks  and  physical  properties. 

The  hygrothermal  conditions  which  were  used  to  condition  the  lam¬ 
inates  were  98  percent  relative  humidity  and  water  immersion  each  at 
room  temperature  and  at  120°,  150°,  180°  and  200°F.  For  each  laminate 
orientation,  comparisons  are  made  between  the  effect  of  each  of  the 
temperatures  and  between  the  two  moisture  condtions.  The  coefficient 
of  diffusion  is  calculated  at  each  temperature  from  the  moisture  de¬ 
sorption  curves  of  saturated  specimens. 


ABSTRACT 


Determination  of  the  Relationship  of  Free  Volume  to  Mechanical  Behavior 
for  an  Epoxy  System  Subjected  to  Various 
Aging  Histories.  (December  1979) 

Darryl  Steven  Ring,  B.S.,  Iowa  State  University 
Chairman  of  Advisory  Committee:  Dr.  W.  L.  Bradley 

An  epoxy  resin  commonly  used  in  advanced  composite  materials  for 
aerospace  applications  was  tested  for  viscoelastic  behavior  after  being 
aged  under  various  cooling  rates  and  steady  state  temperatures.  In 
addition,  the  changes  in  free  volume  for  the  same  aging  conditions  were 
studied.  A  discussion  of  the  correlation  between  free  volume  and 
viscoelastic  response  is  included.  Data  obtained  from  these  experiments 
is  useful  in  contributing  to  an  overall  understanding  of  factors  impor¬ 
tant  to  the  optimization  of  processing  parameters  in  the  manufacture  of 
composite  material  components. 

i  L 

til. 
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ABSTRACT 

Aeroelastic  Tailoring  of  Composite  Materials  (August  1979) 

Jesse  Byron  Rogers,  Jr.,  B.S.,  Mississippi  State  University; 

M.S.,  University  of  Texas 

Chairman  of  Advisory  Committee:  Dr.  Walter  J.  Horn 

The  behavioral  characteristics  of  lifting  surface  structures 
may  be  controlled  through  the  utilization  of  properly  tailored 
composite  laminates.  Variables  such  as  the  number  of  material 
systems  used,  the  number  of  layers,  the  properties  of  the  consti¬ 
tuent  fiber  and  matrix  materials  for  each  system  employed,  and  the 
orientation  and  location  of  each  ply  in  the  laminate  may  be 
considered  in  the  design  process.  The  concept  of  controlling  the 
structural  response  by  optimizing  these  variables  is  called 
aeroelastic  tailoring. 

In  this  report,  the  type  of  structural  response  to  applied  loads 
for  various  composite  laminates  and  the  dependence  of  this  response 
on  the  material  parameters  are  illustrated.  Previous  investigations 
of  aeroelastic  tailoring  are  assessed  and  important  areas  requiring 
further  research  are  identified  and  evaluated. 
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6.  PROFESSIONAL  PERSONNEL  INFORMATION 


6.1  Faculty  Research  Assignments 

Each  participating  faculty  member  is  responsible  for  the  research 
conducted  in  at  least  one  specific  area  of  investigation,  as  shown  below. 
In  addition,  most  serve  as  chairmen  of  one  or  more  of  the  graduate 
advisory  committees  for  M.S.  students  and,  as  such,  direct  their 
students'  research  project.  The  faculty  also  contribute  to  other  research 
activities  on  the  project  by  serving  on  student  advisory  committees, 
through  technical  meetings,  informal  discussions,  and,  in  some  cases, 
through  specific  research  work. 

The  Principal  Investigator  (R.  A.  Schapery)  has  responsibility  for 
overall  technical  direction  and  coordination  and  for  project  management. 

In  addition  he  has  direct  responsibility  for  certain  research  work,  as 
noted  below. 


Faculty  Member/Departmental  Affiliation 

Dr.  Richard  Schapery/Aerospace  and 
Civil  Engineering 

Dr.  Walter  Bradley/Mechanical 
Engineering 

Dr.  Walter  Haisler/ Aerospace 
Engineering 

Dr.  Joe  Ham/Physics 

Mr.  Bob  Harbert/Civil  Engineering 

Dr.  Ken  Jerina/Civil  Engineering 

Dr.  John  Randall/Nuclear  Engineering 
Dr.  Jack  Weitsman/Civil  Engineering 


Primary  Research  Responsibility 

Principal  Investigator  and 
Theoretical  Models  for  Physical 
Aging,  Damage  Growth,  and 
Fracture. 

Physical  Aging  Behavior, 
Delamination  Fracture  Properties. 

Development  of  Finite  Element 
Model s . 

Curing  and  Aging  Studies. 

Duomorph  Gage 

Experimental  Data  Base, 

Mechanical  and  Failure  Property 
Characterization. 

Neutron  Radiography. 

Constitutive  Relations, 
Environmental  Effects. 


6.2  Additional  Professional  Staff 


Mechanics  and  Materials  Research  Center: 

Mr.  Carl  Fredericksen  -  Electronics  Technician 
Mr.  Mark  McEndree  -  Research  Assistant 
Mr.  William  Eue  -  Computer  Programmer 

Nuclear  Science  Center: 

Mr.  Joseph  Taft  -  Research  Engineer 

6.3  Spoken  Papers  and  Lectures  at  Conferences  and  Other  Professional 
Activities  of  the  Faculty  Related  to  Composite  Materials 
(15  January  1979  -  31  December  1979): 


R.  A.  Schapery 


Conference  Presentation: 

"Analysis  of  Damage  Growth  and  Failure  for  Viscoelastic  Composite 
Materials",  Society  of  Rheology  Golden  Jubilee  Meeting,  Boston, 
October  1979. 

Workshops: 

Participant  in  "Workshop  on  Damping  in  Composite  Materials",  MIT. 
Sponsored  by  AFFDL.  April  1979. 

Co-chairman  of  Polymer  Science  Division  of  "Workshop  on  Micro¬ 
structure  Mechanisms  and  Macromechanical  Behavior  of  Materials", 
Cincinnatti,  Ohio.  Sponsored  by  NSF  and  AFOSR.  May  1979. 

Technical  Committee  Membership: 

Member  of  the  committee  "High  Temperature  Metal  and  Cermaic  Matrix 
Composites",  National  Materials  Advisory  Board,  Washington,  D.C. 
Sponsored  by  DOD.  1979-1980. 


W.  L.  Bradley 


Conference  Presentat ions: 

Gordon  conference  lecture  on  research  in  progress  in  the  area  of 
physical  aging  of  epoxy  resin,  Santa  Barbara,  California. 

January  1979. 

"Fracture  Toughness  Studies  of  Cast  Iron  using  a  J-integral  Approach" 
Third  International  Conference  on  Mechanical  Behavior  of  Materials, 
Cambridge,  England.  August  1979. 

Technical  Committee  Membership: 

Member  of  committee  ASTM  E-24  on  fracture. 

Member  of  Metal  Properties  Council,  sponsored  jointly  by  ASME,  ASM 
and  ASTM. 


W.  E.  Haisler 


Conference  Presentation: 

"An  Incremental  Form  of  the  Single-Integral  Nonlinear  Viscoelastic 
Theory  for  Elastic-Plastic-Creep  Finite  Element  Analysis," 

ASME  Pressure  Vessels  and  Piping  Conference,  San  Francisco,  Calif. 
June  1979. 


K.  L.  Jerina 

Technical  Committee  Membership: 

Member  and  Corresponding  Secretary  for  ASTM  Committee  E-9  on 
Fatigue. 

Member  of  the  SAE  Committee  on  Fatigue  Design  and  Evaluation  and 
Chairman  of  the  Task  Force  on  Composite  Materials: 

May  1979:  Organized  and  chaired  a  technical  session  on 

manufacturing  technology  for  automotive  composites. 

Oct.  1979:  Organized  and  chaired  a  technical  session  on 
fatigue  and  design  of  automotive  composites. 

Y.  Weitsman 


Conference  Presentation: 

"Residual  Stresses  due  to  Cool-Down  of  Epoxy  Resin  Composites", 
ASME  Winter  Meeting,  N.Y.  December  1979. 
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In  addition  to  the  above  activities,  the  faculty  attended  several 

i 

conferences  on  composites,  lectured  at  various  Universities  and  Industries 

! 

and  worked  as  consultants  to  industry  in  the  composites  area.  i 

i 
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APPENDIX  A 


Recently  Completed  Reports  and  Publications 


1.  "Interfacial  Stresses  in  Viscoelastic  Adhesive-Layers  Due  to 
Moisture  Sorption"  by  Y.  Weitsman. 

2.  "Residual  Thermal  Stresses  Due  to  Cool-Down  of  Epoxy-Resin 
Composites"  by  Y.  Weitsman. 

3.  "Optimal  Cool-Down  in  Linear  Viscoelasticity"  by  Y.  Weitsman. 

4.  "Residual  Stresses  due  to  Environmental  Conditioning  of  Cross- 
Ply  Graphite/Epoxy  Laminates"  by  Y.  Weitsman. 

5.  "Thermal  Buckling  of  Free  Plates  due  to  Non-Uniform  Temperature 
Distribution"  by  Y.  Weitsman. 

6.  "A  Method  For  Determining  the  Mode  I  Delamination  Fracture 
Toughness  of  Elastic  and  Viscoelastic  Composite  Materials" 
by  D.F.  Devitt,  R.A.  Schapery,  and  W.L.  Bradley. 


hu  J.  Solids  Structures  Vol.  15.  pf>  701-713 
Peryuua  Press  Ltd..  1979.  Printed  ia  Great  Britain 


INTERFACIAL  STRESSES  IN 
VISCOELASTIC  ADHESIVE-LAYERS 
DUE  TO  MOISTURE  SORPTION 

Y.  Weitsman 

Mechanics  and  Materials  Research  Center,  Texas  A&M  University,  College  Station,  TX  77843,  U.S.A. 

( Received  30  August  1978;  revised  IS  December  1978) 

Abstract— This  paper  concerns  the  mechanical  behavior  of  an  epoxy  adhesive  layer  that  is  located 
between  stiff  adherends  as  the  adhesive  absorbs  moisture  from  the  ambient  environment. 

As  the  moisture  penetrates  the  layer  by  an  extremely  slow  diffusion  process  the  eposy  undergoes  a 
simultaneous  process  of  stress-relaxation.  Calculations  of  viscoelastic  interfacial  stresses  were  performed 
for  a  time-dependent  response  which  is  typical  of  epoxy  for  a  layer  geometry  as  encountered  in  practice. 

The  results  show  that  for  exposure  to  steady  ambient  humidity  the  viscoelastic  stresses  are  smaller  than 
their  elastic  counterparts.  However,  under  fluctuating  ambient  humidity  the  viscoelastic  response  may 
cause  stress  reversals,  and  thereby  failure  modes,  which  are  not  predicted  by  elasticity  theory. 


1.  INTRODUCTION 

Adhesive  bonding  forms  an  attractive  method  of  joining  structural  members.  In  analyzing  the 
stresses  which  develop  within  the  bond  it  is  necessary  to  account  for  the  fact  that  the  adhesive 
materials  respond  in  a  viscoelastic  manner  under  loads  and  their  time-dependent  behavior  is 
greatly  affected  by  moisture  and  temperature. 

When  two  adherends  are  joined  together  by  a  thin  adhesive  layer  the  adhesive  absorbs 
moisture  from  the  ambient  environment,  at  its  exposed  edges,  which  induces  swelling  strains 
into  the  layer.  Since  the  adherends  are  much  stiffer  than  the  adhesive  they  constrain  the 
adhesive-layer  against  its  tendency  to  swell,  thus  causing  the  formation  of  extremely  high 
stresses  within  the  layer.  In  this  paper  attention  is  focused  on  the  interfacial-stresses  which 
arise  at  the  interfaces  between  the  adherends  and  the  adhesive. 

Since  the  moisture  penetrates  the  layer  by  an  extremely  slow  diffusion  process,  the  epoxy 
may  undergo  noticeable  creep  and  relaxation  while  the  diffusion  process  is  in  progress.  The 
main  purpose  of  this  paper  is  to  relate  the  interaction,  which  occurs  concurrently,  between  the 
two  time-dependent  phenomena — moisture-diffusion  and  stress-relaxation. 

The  analysis  employs  a  variational  method  and  is  inherently  approximate  in  nature.  It  is 
due  to  this  approximation  that  the  edge  singularity  in  the  stress  field  is  replaced  here  by  finite, 
though  large,  values.  However,  for  the  exceedingly  thin  layers  that  are  utilized  in  practice,  the 
selected  expressions  for  the  displacement  fields  should  provide  a  good  approximation.  Fur¬ 
thermore,  the  approximation  should  not  detract  from  the  basic  purpose  of  this  work  which  is  to 
provide  information  about  the  relative  influence  of  the  diffusion  and  relaxation  times. 

2.  FORMULATION 

Consider  a  semi-infinite,  isotropic  and  elastic  adhesive  layer,  of  thickness  2a,  between  two 
semi-infinite  adherend  plates  as  shown  in  Fig.  1. 

Let  X,  Y  denote  Cartesian  coordinates  and  t  time.  Let  e(X,  t )  represent  the  unconstrained 
swelling  induced  by  moisture  and  n,  v  be  the  elastic-shear-modulus  and  Poisson's  ratio, 
respectively.  Assume  a  state  of  plane-strain. 

Introduce  the  following  non-dimensional  quantities 

x  -  Xla,  y  =  Yla,  e(x,  t)  =  e(X/a,  t) 

.  .  (1) 
s,  =  s,  =  a,  fi,  s„  =  t„Iix. 


Also,  denote 


/,  =  2W(I  -  2v),  U  =  2(1  -  v)/(l  -  2v),  /j  =  2(1  +  „)/(!  -  2v). 
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Rg.  1.  The  adhesive  layer  between  rigid  adhere nds  exposed  to  moisture  at  X  =  0. 


Consider  now  the  following,  approximate,  form  for  the  displacement  fields  within  the 
adhesive-layer 


U(X,Y)  =  au*x)+\au4x)yi 

,  (2) 

V(X,  Y)  =  a  v,(x)y  +ja  o, (x)y1. 

In  (2)  U  and  V  are  components  of  displacement  in  the  X  and  Y  directions,  respectively, 
and  Uo,  tt2,  ft,  t>3  are  dimensionless,  yet  to  be  determined  functions  of  the  dimensionless 
coordinate  x.  Note  that  for  an  induced  swelling-strain  of  the  form  e{x,  t)  the  displacements  U 
and  V  possess  the  required  symmetries  in  Y. 

Form  (2)  is  the  lowest-order  approximation  which  provides  information  on  the  interfacial 
stresses,  at  the  boundary  between  adhesive  and  adherend,  in  a  self-contained  manner.  For 
exceedingly  thin  layers  the  approximation  should  provide  sufficient  accuracy. 

Expressions  (2)  are  essentially  the  same  as  used  in  Ref.  [1]. 

Employing  Hooke's  law  we  obtain  the  following  non-dimensional  form  for  stress-strain 
relations: 


S* 

'h  /.  U  0 

Sf 

= 

/.  h  U  o 

Sir 

0  0  0  1 

_  - 

ut'  +  ^ut'y1 

v,  +  v,y2 
-e 

u2y  +  v  ,'y  +  jv}'y3 


(3) 


where  primes  designate,derivatives  with  respect  to  x. 

In  addition,  due  to  the  assumed  rigidity  and  immobility  of  the  adherends  we  have  the 
following  boundary  conditions: 


n(x,  l,f)*0 
p(x,  1,0  =  0. 


(4) 


In  view  of  the  constraints  implied  by  boundary  conditions  (4)  it  is  obvious  that,  when 
employing  a  variational  formulation,  the  unknown  functions  a*  u2,  t>i.  and  Pj  in  (2)  cannot  be 
varied  independently.  As  noted  in  [1]  the  four  above  mentioned  functions  can  be  treated 
independently  of  each  other  by  introducing  suitable  Lagrange  multipliers  >,  and  rh  It  turns  out 
that  s,  and  s„  which  are  the  dimensionless  force-conjugates  of  p(x,  1,0  and  u(x,  1,0  respec¬ 
tively,  represent  the  sought  after  interfacial  normal  and  shear  stresses. 

Employment  of  the  principle  of  virtual  work,  and  utilization  of  the  Lagrange  multipliers  in 
analogy  with  [l|,t  yields  the  following  expression  for  the  variation  of  the  internal  energy  SE 


tSee  particularly  the  aectkw  Variational  Formulation  of  the  Totally  Constrained  Cue  (u»-0). 
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SE  =  tia2  |j^  1  y28u}')  +  s,(6v,  +  y*Sc })  +  s„(y««2  +  y«e,'  + j  yJ«»j')J  dy 

-  s, ( 6ii<,  +  ^  +  |  u2j  -  sM(sv,  +j  8v}j  -  8s.(v,  +  j  dx.  (5) 

Similarly,  the  variation  of  the  external  work,  8W,  is  given  by 

8W  =  tia2  /#[ix*(0,  y)(&io  +  \  y’Su  2)  +  r„*(0,  y)(y8v,  +  ~  y3St>3)]  dy.  (6) 

Integration  by  parts  of  (5)  and  employment  of  (3)  yield,  upon  factoring  out  the  now 
independent  variations  8u o,  8u2,  8v,  and  Sv2: 

\fitf-Ue'  +  s,  =  Q 

JS  htf  ~  fat  +  l)o,'  + 1  «o  - 1  he'  +  \  s,  =  0 
2  2  (7) 

f(/s  +  l)«o'  -  ^  0,"  - 1  /20 1  -  he  -  i«  =  o. 

To  these  are  adjoined  the  constraint  conditions 

«o  +  ^«2  =  0 

(8) 

v,  +103  =  0. 

Equations  (7)  and  (8)  express  the  field  equations  of  the  problem.  The  concomitant  from  the 
integration-by-parts  of  (S)  combines  with  (6)  to  determine  the  boundary  conditions  at  x  =  0. 
This  combination,  together  with  (3),  yields 


-7«o  +  2o,'  =  0 

2«o'  +  ^oi  =|^e(0,  t). 


(9) 


Note  that  the  boundary  conditions  (9)  were  obtained  after  employment  of  (8),  which  hold 
also  at  x  -  0. 

In  order  to  solve  for  u0  and  t>,  we  first  eliminate  s,  between  the  first  pair  of  (7)  and  s, 
between  the  last  pair.  Next,  we  employ  (8)  to  eliminate  u2  and  p3.  This  procedure  yields: 


2h«o'+(h  +  l)vl'-S*o  =  lhe' 
j  ®r~(/i  +  1)po' ~3/2p,  =  0. 


(10) 


The  solution  of  (10)  is  expressed  with  the  aid  of  the  following  second  order  equation 


jh?  +  [(/.  +  «*  ■ ~  6/a*  -  ^]z  +  15  h  =  0. 


(11) 


I 
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Denote  the  roots  of  (II)  by  Zi  and  Zj  and  define 


—  Vzi,  ati  -  'fZi 


-a? 


x  =J£i±Jl2y_  D=ifi - ! - 

*U  (2/7)Zu-3 h  *h»iXi-*Xxm 

The  solution  of  (10),  which  vanishes  as  x -*<*>,  is 

u0=  D  j  sinh  aj(x  -  r)«'(r) ds~xtj  sinh  a,(x  -  s)e'(s) dr  j  +  At  e'*1*  +  A2 e' 

vt  =  DxiXif  [cosh  ajx  -  r)  -  cosh  a,(x  -  r) j  e'(s)ds  -  XiA,  e~m,‘ -  x2 A2e~°11‘.  (12) 

In  (12)  A,  and  A2  are  arbitrary  constants  that  are  determined  with  the  aid  of  boundary 
conditions  (9). 


3.  ELASTIC  SOLUTION  FOR  MOISTURE  DIFFUSION  UNDER  CONSTANT  AMBIENT  HUMIDITY 

It  has  been  observed  [2],  [3]t  that  the  saturation-moisture  level  in  “neat”  epoxy,  as  well  as  in 
epoxy-based  composites,  depends  on  the  relative  humidity  of  the  ambient  environment.  The 
dilatational  strains  that  accompany  moisture  sorption  are  about  2/3  of  the  swelling  that  would 
be  anticipated  by  straightforward  “volume  additivity”. 

The  penetration  of  moistures  into  several  epoxy  resins  was  shown  to  follow  the  classical 
diffusion  equations  (2-4]. 

For  the  one-dimensional  diffusion  process  considered  herein,  it  follows  that  under  constant 
ambient  humidity  the  moisture-induced  swelling  is  given  by 


e(X,O  =  /terfc0=).  (13) 

In  (13)  A  is  a  constant  which  converts  moisture  content  to  dilatational  strain  and  k  is  the 
coefficient  of  moisture  diffusion. 

Introduce  the  non-dimensional  time  (3  =  to/a2,  then  the  non-dimensional  version  of  (13) 
reads 


e(x,  P)  =  A  erfc(x/2V^).  (14) 

Substitution  of  (14)  into  (12)  leads  to  a  closed-form  solution  for  all  the  integrals  therein  [$,  6], 


eg- 


J  sinh  a(x  -  r)e'(r)  dr  =  “e^’[sinhox[l-^erf^^=  +  aV0^ 

- "im ''*))]■ -  2  c“h  + * -  Aw ■  °  ^)] } 


and  a  similar  expression  for 


j  cosh  o(x  -  s)e'(s)  dr. 


With  the  simpler,  closed-form  version  of  h«  and  t>(  that  replaces  all  the  integrals  in  (12)  it  is 
possible  to  evaluate  analytically  the  unknowns  At  and  A2  therein. 


tFor  most  of  (Ms  information,  the  author  it  (ready  indebted  to  Metart.  J.  E.  Halkias  and  E.  L.  McKafue  of  General 
Dynamics  Corporation,  Fort  Worth  Division. 
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Employing  the  boundary  conditions  (9),  with  e(0,  i)  =  A,  the  values  of  A  i  and  At  are  determined 
from  the  following  two-by-two  algebraic  system: 


L,  L{ 

M 

M,  Mi 

i 

Ah 


D(~  L-iXiJi  +  L[Xih) 
A  +  D(-  M\XiKt  + 


(15) 


where,  in  (15) 


Lu  =  7-(2ax)a  Mu=  (^a  +^)u 

Ku  =  -Ae'k,’J  A.2=-Kuerf(auV0). 

Solving  for  At  and  A2  from  (15)  we  obtain  the  complete  solution  to  u0  and  t>,  in  (12).  This,  in 
turn,  enables  us  to  evaluate  the  interfacial  tractions  s„  and  s,  by  using  the  first  and  third  of  eqns 
(7). 

After  straightforward,  though  tedious,  manipulations  we  obtain 


s„(x,  0)  =  |d  [a2X\P2KJ.x,  P)  -  cttxi  PiKt(x,  P)] 

- 1  a,P,  A,  e'”'1  -  \  a jP2A2  t  -  /,  e(x,  p) 

s,(jc,  P)  =  - 1  ftD  latxM**  P)  -  afxJifa  P)  1 

-|/Wi*,e-^  +  4h,%e-^)-j^e-««. 


(16) 


In  (16) 

/.i(*.«--Y^,,u(#w-#Tj) 

where 

d'l  j  =  e'^erfctauV^  +  xl2Vp) 
tii  =  e-'^erfc  (~auVp  +  xllVp). 

Also, 

P\2  =  fl  +  l~|(o*,)u-  i 

Form  (16)  is  the  elastic  solution  to  the  present  problem.  While  this  form  provides  a 
necessary  step  toward  the  generation  of  the  viscoelastic  solution,  because  its  Laplace  transform 
is  readily  available,  it  cannot  be  evaluated  numerically  for  a  wide  range  of  x  and  t  because 
Ji(x,p),  K,(x,  p),  A i  and  A2  become  numerically  unstable.t  Difficulties  of  this  sort  stem  from 
the  peculiar  geometry  of  the  adhesive  layer,  namely,  the  extremely  small  ratio  of  thickness  to 
length. 


tPor  x  -  3  and  t  -  1000,  with  typical  values  for  k,  »  and  »,  /,  and  K,  arrive  products  of  Bombers  of  0(10“)  with 
numbers  of  0(10'”).  thus  esceedmt  the  accuracy  provided  even  by  "double  precision”. 
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The  difficulty  is  overcome  by  the  employment  of  an  approximate  form  of  erfc  (z),  as 
follows  [7]: 

erfc  2  =  Q(t(z))e-‘l 

where  (l7) 

s 

QU)sj£  and  f(r)  =  (l+pzr'$. 

Appropriate  substitution  of  (17)  into  (16)  yields  an  alternative,  numerically  stable,  form  for 
the  interfacial  stresses.  We  obtain 

J  s„(x,  fi)  - 1  ^  -  u  tdcmVfi) 

-  H(x  - latf) e^-*»]  +  «,*2P|  [j  e  ,xl,4',(<3(a,V^  + 

+  h(x  -  2a,p)Q  (|  a,Vfi-  ^|))  -  H(x  -  2a  ip)  e'*-'2-'] 

+  |(a,/W2  e"'*  -  a2P2Mt  e"*)} 


i  „U,K)  -  a, ■L,.--) < 

*  *  3j)  -  *<*  -  **»(  |«.V?  -  ^jl)) 

+  W(jr  -  la-if})  e**-**]  -  e^J,w(<?(a,V^  +  ^=) 

-  A(*  -  2a,0)o(|a,Vi3  -  ^=|))  +  H(x  -  2a, 0)  e'ta'J"’*] 

+  f(a,iM1e-«'I-«JlM,e  •>*)}. 


In  (18) 


«<»-{?  fill  *«>-{'!  !>! 


*  =  LxMi~  LiM\,  B  - X2L]QialVp)-x,L2Q(a2Vp). 

As  f-»®  the  “fully  saturated”  elastic  solution,  obtained  by  taking  limits  of  (18)  as  £-♦«>, 
reads 

»-(*.*)  =  ~/j  f1  + 1  ^(a,.P,/.2  e -|X-  a2P2L\  e'**)] 

1  5/  (,9) 

^  •»<(*. »)  =  -  j  ^  (a,2L2  e  -  a2lZ„  e  -•**). 

file  eemencei  value  of  a,  (/«  1,2,.  ..5) and p  are  given  in  Section  7. 1 .26,  on  page  299  of  Ref.  [7]-  This  approximation  is 
valid  for  all  r. 
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4.  VISCOELASTIC  SOLUTION  FOR  MOISTURE  DIFFUSION 
UNDER  CONSTANT  AMBIENT  HUMIDITY 

Uniaxial-tension  experiments  {?]  on  epoxy  indicate  that  its  creep-response  can  be  described 
by  a  “power  law”  form  as  follows: 

D(l;r)=A,(r)  +  A[^yj]*  (20) 

In  (20)  T  denotes  temperature  and  a(T)  is  the  “shift-factor”  function.  Additional  data  on 
polyester  [9]  and  epoxy  [10]  indicate  that  the  creep-response  depends  also  on  the  moisture 
content  m,  namely  D  =  D(t;T,m).  However,  this  dependence  is  omitted  in  the  present 
analysis.  Furthermore,  we  shall  consider  the  case  of  uniform  temperature  for  the  entire 
adhesive. 

In  most  adhesives  the  Poisson’s  ratio  v  remains  constant  under  a  wide  range  of  conditions. 
Consequently  we  shall  consider  a  viscoelastic  shear-compliance  of  the  form 

p(0  =  d„(l  +  dt/«) 

where 

do  =  2(1  +  v)  A>.t  (21) 

Denote  r  =  a2 Ik  and  r,  =  df ,/*.  Let  p  be  the  variable  of  the  Laplace  transform  in  the  time 
domain  and  ij  =  pr.  Finally,  let  the  ratio  between  the  characteristic  diffusion  time  r  and  the 
characteristic  creep  time  t,  be  6  =  tJt. 

Then,  the  Laplace-transform  of  (21)  is 

The  reciprocity  relations  between  the  transforms  of  the  creep  and  relaxation  functions[U] 
yield 

d<>pE(p)  =  l/[  1  +  T(q  +  »/(*»)«]  =  R«h >).  (23) 

Upon  insertion  of  the  explicit  expressions  for  A,  and  A2  into  (16),  the  Laplace  transforms  of 
the  elastic  stresses  (denoted  by  f )  become 


/  auftPi  e~v™-e-°'11  a-iX\Pi  e^-e^X 
'Vij+ai  Vrj-ctt  Vij+a  2  y/rj-a2  / 
XiLi  X\L2  \  a,P,M2e~",x  -  a2P2M\  c-***1'] 
'Vi)  +  ai  Vrj+ay  A  J 

-vs.  ■  5/iatFiL,e-^-«,P,L,e~a|J 


ftLi  _  X\L2  \  atiP,M2e~",x  -  a2P2M,  e'^l 
Tf  +  at  \/ti+a  -J  A  J 

.  5  U  a2PiLi  e-***  -  a,P,L2  c'”'1 
3  h  A 

-vj, _ 5  - 

3J  A 

-r  «i  lXi  Vrie^'-aie-^*  g22x,  VitC^-gyC"^* 
LVij+ai  y/rj-at]  y/rj+a2  y/rj-a2 

ri  y/rj  +  ay  A  -I 

rities  at  Vi j  =  a,  and  Vi)  =  a2  are  removable, 
elastic  stresses  are  liven  by  s„  » <r,Jd<y 
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From  the  fundamental  equations  of  linear  viscoelasticity  [12]  it  follows  that  the  Laplace 
transforms  of  the  viscoelastic  solution  to  our  problem  (denoted  by  overbars)  are 

j  psn(x,  p)  =  R(Bt})  •  ^  psm(.  x,  p) 

,  .  <25> 

^ps,(x,  p)  =  R(0ri)  ■  —  ps,(x,  p ). 

In  order  to  obtain  the  viscoelastic  solution  in  the  real-time  domain  it  is  necessary  to  invert 
(25).  The  inversion  was  obtained  by  the  method  of  collocation  [13],  Typically,  expressions  (25) 
were  plotted  vs  log  r\  and  a  set  on  N  points  ij,(i  =  I, . . .  N)  was  selected  to  cover  the  ranges 
where  those  plots  show  noticeable  curvatures.  It  was  found  that  the  shapes  of  the  above 
mentioned  plots,  as  well  as  the  locations  of  high  curvature  regions,  depended  on  x.  Therefore,  it 
was  impossible  to  employ  any  specific  set  of  points  i?(  in  the  inversion  and  an  appropriate 
selection  had  to  be  made  for  each  value  of  x. 

The  calculations  followed  the  scheme  of  Ref. [1 3).  Note,  however,  that  for  a  “power-law” 
creep  as  given  in  (20)  we  have  lim  D(t,  T)  -*  <*>,  hence  the  equilibrium  stress  vanishes.  Further- 

more,  in  the  present  case  the  transform-variable  is  r>  =  pr,  whereby  the  inversion  is  obtained  in 
terms  of  dimensionless  time  th,  rather  than  real  time  /. 

Calculations  were  performed  for  an  adhesive-layer  with  the  following  selected  properties: 

Thickness  a  =0.003",  moisture  diffusivity  k  =  0.2  x  10~8  in2/sec,  Poisson’s  ratio  v  =  0.45, 
100%  R.  H. -to-swelling  conversion-factor  A  -  0.03,  creep  compliances  D0  =  2  x  1(T6  in2/lb, 
Di  -  0.2  x  I0-6  in2/lb  with  q  -  0.19.  The  corresponding  value  of  the  shift-factor  function  was 
a(T)=  1. 

With  the  above  numerical  values  we  obtain  t  =  4500  sec-1  and  the  value  of  6  in  (22)  is 
8  =  2450. 

Numerical  evaluations  were  performed  on  an  Amdahl  470  digital  computer.  The  results  are 
shown  in  Figs.  2-9.  The  heavy  solid  lines  in  Figs.  2-7  show  the  dependence  of  the  non- 
dimensional  viscoelastic  normal  and  tangential  interlaminar  stresses  s„  and  s,  on  time  t  at  three 
fixed  stations  within  the  layer,  x  =  XI a  =  0.01, 1  and  3,  when  the  free  edge  (x  =  0)  is  exposed  to 
a  constant  R.H.  =  100%  at  t  -  0.  Note  that  the  time  scale  is  logarithmic. 

For  comparison  purposes,  the  elastic  stresses — as  given  in  (18) — are  shown  in  thin  solid 
lines  in  Figs.  2-7.  Note  that  any  station  within  the  layer  senses  the  build-up  and  approach  of  the 


Fig.  2.  Elastic  and  viscoelastic  values  of  the  non-dimensional  normal  interlaminar  traction  »,  at  x  =  XI a  - 
0.01  vs.  log  /  (I  in  seconds).  Heavy  lines— viscoelastic,  thin  lines— elastic.  Solid  lines— exposure  to  a 
constant  ambient  R.H.,  dashed  lines— exposure  to  fluctuating  ambient  R.H.  (according  to  eqn  29). 
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Fig.  3.  Elastic  and  viscoelastic  values  of  the  non-dimensional  tangential  interlaminar  traction  s,  at 
x  =  Xla  =  0.01  vs.  log  /  (1  in  seconds).  Heavy  lines— viscoelastic,  thin  lines — elastic.  Solid  lines— exposure 
to  a  constant  ambient  R.H.,  dashed  lines— exposure  to  fluctuating  ambient  R.H.  (according  to  eqn  29). 


Fig.  4.  Elastic  and  viscoelastic  values  of  the  non-dimensional  normal  interlaminar  traction  s.  at  x  «  Xla  =  I 
vs  log  t  (t  in  seconds).  Heavy  lines— viscoelastic,  thin  lines— elastic.  Solid  lines— exposure  to  a  constant 
ambient  R.H.,  dashed  lines— exposure  to  fluctuating  ambient  R.H.  (according  to  eqn  29). 


moisture  in  a  gradual  manner  and  that  a  steep  increase  in  stress  occurs  after  the  passage  of  a 
time-span,  which  depends  on  the  characteristic  diffusion  time  and  the  location  of  the  station. 
Beyond  that  time  the  station  notices  a  “fully  saturated”  state,  which  explains  the  leveling-off  of 
the  elastic  stresses  for  longer  times. 

The  viscoelastic  values,  which  are  affected  by  relaxation,  diminish  to  zero  with  time. 

Thus  far  all  viscoelastic  results  were  obtained  from  inversion  of  (25)  by  means  of  collo¬ 
cation. 
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fig.  S.  Elastic  and  viscoelastic  values  of  tie  non-dimensional  tangential  interlaminar  traction  s,  at  Jt  =  XI  a  -  I 
vs  log  /  «  in  seconds).  Heavy  lines— viscoelastic,  thin  lines— elastic.  Solid  lines— exposure  to  a  constant 
ambient  R.H.,  dashed  lines— exposure  to  fluctuating  ambient  R.H.  (according  eo  eqn  29). 


Fig.  6.  Elastic  and  viscoelastic  values  of  the  non-dimensional  normal  interlaminar  traction  i„  at  x  =  XI a  =  3  vs 
log  l  (r  in  seconds).  Heavy  lines— viscoelastic,  thin  lines— elastic.  Solid  lines— exposure  to  a  constant  ambient 
R.H.,  dashed  lines— exposure  to  fluctuating  ambient  R.H.  (according  to  eqn  29). 


Consider  now  an  approximate  expression  for  the  relaxation  modulus(14] 


where 


E(t) 


1  sin  nn 
D(t )  m 


(26) 


riSf?0*"- 

Employing  (26),  the  “quasi-elastic”  (viscoelastic)  solution  is  given  by(13] 

*?*■(*,  0  ■  F(tUm(x,  t) 

*?*(*,  0  -  FV)s,(x,  f) 
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Fig.  9.  Elastic  and  viscoelastic  values  of  the  non-dimensional  tangential  interlaminar  traction  s,  at  times 
I  =  1000  sec  and  f  -  10,000  sec  vs  x.  Dashed  lines— elastic,  solid  lines— viscoelastic. 

Employing  (27),  values  of  and  s‘‘(x,  t)  were  calculated  for  fixed  values  of  time, 

(  =  1000  sec  and  10,000  sec,  along  the  interface  between  the  adhesive  and  adherend.  These 
results  are  shown  in  Figs.  8  and  9  for  0  s  x  s  1.6.  (For  comparison  purposes,  the  elastic  values 
are  also  shown  in  dashed  lines.) 

5.  VISCOELASTIC  SOLUTION  UNDER  FLUCTUATING  AMBIENT  HUMIDITY 

Thus  far  the  adhesive-layer  was  considered  to  be  subjected  to  a  sudden  exposure  to  ambient 
humidity  at  time  t  =  0,  with  a  swelling  conversion  factor  A. 

For  fluctuating  ambient  humidity  A  =  A(t)  the  interlaminar  stresses  are  given  by  a  super¬ 
position  integral  as  follows 


s'S(x,t-n 


dA(Q 

dt' 


dt' 


(28) 


where  m  =  n  or  /. 

In  (28)  the  kernel  functions,  under  the  integral  sign,  are  given  by  (27). 

A  sample  computation  was  performed  for  a  sudden  exposure  to  100%  R.H.  followed  by 
drying  at  time  t  =  lt  later.  In  this  case 

A(f)  =  A[tf(f)-  H  (/-/,)]  (29) 


whereby 


»«*(*,  t)  *  s’, ;*fx,  t) -  s’„(x,  t-t,)  (m  =  n  or  /).  (30) 

Results  for  t,  =  10,000  sec.  at  x  =  0.01,  x  =  1.0  and  x  =  3.0  are  shown  by  the  dashed  fines  in 
Figs.  2-7.  Note  that  the  superposition  of  wetting  and  subsequent  drying  is  algebraically  additive 
in  the  elastic  case.  Therefore,  all  the  elastic  stresses,  represented  by  thin  dashed  lines  in  Figs. 
2-7,  tend  to  zero  as  /  -* »  without  any  further  sign  changes. t  However,  in  the  viscoelastic  case, 


tThere  it  nevertheless  in  exception  m  Fig.  3. 
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the  drying  effects  are  superimposed  on  a  relaxed  wet  state  and  this  causes  an  overcompen¬ 
sation— thereby  introducing  sign-reversals  in  the  stresses.  Eventually,  as  f-»«,  the  viscoelastic 
stresses  too  tend  to  zero. 

It  is  observed  that  viscoelasticity  predicts  detrimental  effects  that  are  caused  by  fluctuations 
in  relative  humidity— such  as  low  cycle  fatigue  and  interfacial  tension— that  are  absent  in  an 
elastic  analysis. 


6.  CONCLUDING  REMARKS 

The  present  analysis  is  based  on  several  idealizations  which  restrict  the  validity  of  the 
results.  The  exceedingly  large  stresses  which  are  predicted  by  linear  theories  near  the  corners 
X  =  0,  Y  =  ±a  cannot  be  borne  by  real  adhesives  and  the  response  in  those  regions  must  be 
represented  by  a  suitable  non-linear  model. 

In  addition,  the  actual  viscoelastic  response  of  resins  depends  on  the  moisture  content, 
which  introduces  yet  another  non-linearity  into  the  analysis  of  adhesive  layers. 

The  evaluation  of  those  non-linear  effects  still  awaits  an  appropriate  characterization  at  the 
present  time.  Once  available  it  would  require  solution  by  means  of  numerical  methods. 
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Residual  Thermal  Stresses  Due  to 
Cool-Down  of  Epoxy-Resin 
Composites 


This  paper  concerns  the  residual  thermal  stresses  that  develop  within  the  resin  of  a  fiber  - 
reinforced  composite-laminate  as  the  material  is  cooled  from  cure  temperature  down  to 
room  temperature.  The  calculations  presented  herein  consider  the  viscoelastic  response 
of  the  resin  and  account  for  the  temperature  dependence  and  the  stress  sensivitity  of  the 
creep  compliance.  Comparisons  with  linear  elasticity  indicate  that  viscoelastic  relaxation 
may  reduce  the  residual  stresses  by  about  20 percent. 


Introduction 

Linear  elastic  stress  analysis  predicts  that  the  thermal  stresses, 
which  develop  within  a  composite  material  due  to  typical  temperature 
excursions  experienced  by  the  structure,  may  exhaust  the  strength 
of  various  laminae  even  before  .the  application  of  external  loads.  Ef¬ 
fects  which  were  attributed  to  high  residual  thermal  stresses  in 
composite  laminates  were  reported  by  several  investigators  [l,  2j.  A 
similar  prediction  can  be  made  for  the  micromechanical  thermal 
stresses  within  the  epoxy  resin  when  its  response  is  assumed  to  follow 
linear  thermoelasticity.  The  assumption  of  linear-elastic  behavior 
tends  to  overestimate  the  residual  thermal  stresses  and,  by  ignoring 
the  time,  temperature,  and  stress  dependencies  of  the  resin’s  response 
rules  out  any  advantage  that  may  accrue  from  the  consideration  of 
special  temperature  histories. 

A  major  portion  of  the  thermal  stress  is  introduced  into  the  com¬ 
posite  laminate  during  cooling  from  cure  temperature  down  to  room 
temperature.  This  cool-down  stage  may  be  viewed  as  part  of  the 
manufacturing  process  and  therefore  lends  itself  to  careful  control. 
It  has  already  been  shown  that  the  selection  of  an  optimal  cool-down 
path  could  reduce  residual  thermal  stresses  in  plexiglass  |3)  and  the 
same  idea  can  be  carried  over  to  epoxy  and  epoxy-based  compos¬ 
ites. 

This  paper  presents  calculations  of  the  "micromechanical”  residual 
stresses  within  the  epoxy  due  to  cool-down.  Our  purpose  is  to  provide 
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a  quantitative  assessment  of  the  role  of  inelastic  response.  The  cal¬ 
culations  are  based  upon  a  therroorheological  ly  complex  viscoelastic 
model,  which  includes  some  aspects  of  the  nonlinearity  that  is  at¬ 
tributed  to  stress  effects  on  relaxation. 

The  findings  of  this  paper  are  preliminary  since  they  are  based  upon 
incomplete  data.  Consequently  it  was  necessary  to  extrapolate  existing 
data  to  higher  temperature  ranges  and  to  infer  the  nonlinear  stress- 
effects  by  analogy  from  similar  materials  (glass-reinforeced  epoxy). 
Perhaps  most  significantly,  due  to  the  complete  lack  of  data,  it  was 
necessary  to  omit  any  aging  effects,  which  may  be  very  prominent 
during  the  cool-down  of  a  newly  created  laminate. 

Several  experimental  programs,  currently  in  progress,  should 
provide  information  which  would  enable  us  to  reassess  and  refine  the 
model  which  is  used  in  the  present  work. 

Formulation 

Consider  a  fiber-reinforced,  multidirectional,  composite  laminate 
-  h  <  z  <  h,  —  °°<  x,y  <  <■>,  traction-free  at  the  boundaries  z  -  ±h 
and  subjected  to  spatially  uniform  temperature  fluctuations.  For 
sufficiently  slow  fluctuation  of  the  ambient  temperature  and  suffi¬ 
ciently  thin  layers  it  is  possible  to  discard  the  transient  states  and 
consider  the  quasi-static  case  of  uniform  temperature  9(z,t )  s  fH±h,t) 
=  #(().' 

Since  the  present  analysis  aims  at  investigating  the  thermo- 
viscoelastic  response  of  the  resin  let  us  assume  that  the  multidirec¬ 
tional  fibers,  all  of  which  lie  in  planes  z  ■  const.,  affect  the  behavior 
of  the  resin  through  the  following  geometric  constraint 

c,  “  ( v  =  cx/A9  (1 ) 

In  ( 1 1  Of  is  the  coefficient  of  thermal  expansion  of  the  fibers,  \Q  is 


1  In  typical  laminates,  the  transient  effects  may  be  ignored  up  to  temperature 
fluctuations  of  5°K/min. 
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the  difference  between  the  initial  and  final  temperatures,  AO  =  0U  — 
8f,  and  <„  ty  are  the  in-plane  strains. 

If  we  assume  in-plane  isotropy  of  the  layer,  and  consider  linear- 
elastic  response,  we  obtain  the  following  results  [3]: 

a,  =  oy  =  <7  =  -cGA8 

(2) 

at  =  0  and  all  shears  vanish. 

where 

1  +  v 

c  =  2 - («„-«/).  (2') 

1  -  f 

In  (2)  and  (2')  G, »  are  the  shear  modulus  and  Poisson's  ratio  of  the 
resin,  and  am  is  its  coefficient  of  thermal  expansion. 

If  the  resin  responds  viscoelastically  equation  (2)  takes  the  form 

<r(f/)«-c  f''Gtf<(/)-f(f')|0<f')df'  (3) 

0 


where  £  is  a  “reduced  time,"  given  by 


ds 

arM 


(4) 


and  C(t)  is  the  shear-relaxation  function  for  the  resin. 

The  purpose  of  the  subsequent  investigation  is  to  consider  and 
evaluate  the  significance  of  various  factors  which  affect  the  function 
«r.  and  investigate  the  influence  of  the  history  0(f)  on  ait/). 


A  Linear  Thermorheologically  Complex  Model  for 
Epoxy 

On  the  basis  of  uniaxial  tension  data  obtained  for  “shell  58-68R” 
epoxy  resin  |4],  for  temperatures  ranging  between  20° F  and  160°F, 
it  is  possible  to  express  the  creep  compliance  as  follows: 


D(6,t)  =  U0(8)  +  D, 


'  t 


(5) 


where 


Do(0)  =  [0.2731  +  0.000574  (0  -  297)]  10~« 

16480  \ 

O,  =  0.01  X  10-«,  «r(0)  =  exp  I—  -  21.821,  q  =  0.19 


In  (5),  D  is  in  (kPa)-1,  the  time  t  is  in  minutes,  and  the  temperature 
0  is  in  degrees  Kelvin.  Furthermore,  the  Poisson  ratio  for  epoxy  is 
approximately  constant  v  »  0.35.  For  temperatures  below  glass 
transition  the  coefficients  of  thermal  expansion  yield  that 


a  =  a„  —  «/  =  5  X  10  5m/m/l°K. 


Consider  a  composite  cooled-down  between  a  cure-temperature 
0o  “  438°K  and  room  temperature  8/  *  298°K.  For  lack  of  more 
complete  data  assume  that  (5)  can  be  extrapolated  up  to  the  elevated 
cure-temperature.  In  terms  of  the  imposed  thermal  strain  («m  - 
«/)(0 o  -  0(f))  and  the  creep  compliance  (5)  we  obtain  the  following 
integral  equation  for  the  stress  a: 

-JL-l8o-0(t)]  =  Do(0)<T(t)  +  Dl  f‘  [{(f)  -  ^n]qHt')dt'  (6) 
l-i/  Jo 


where  «  =  am  -  a/. 

If  the  temperature  drops  linearly  between  0O  and  0/  over  a  time 
interval  </  then 


0(f)  80  +  Rt  with  R  =  (0/  -  0O )/(/ 

and  equation  (6)  reads 
a 


1  -  v 


Rt  *  Do(0o  +  Rt)o(t) 


+Dl  r*[  —  r 

Jo  l Jo  ar(0(.i ))  Jo 


_ds_l" 


<t  (t’)dt'  (7) 


where  in  (7) 


I  6480  \ 

a^(8,)*MPfe-21-82) 


Fig.  1  Dependence  ot  the  residual  stress  a  on  cooling  time  I,;  linear  ther- 
morheoiogically  complex  model.  Linear  (constant-rate)  cool-down  tram  8, 
=  438  K  to  8 1  =  298°K.  Dashed  Nnes  Indicate  further  relaxation  of  residual 
stresses  at  0,  =  298°K  alter  termination  ol  cool-down. 


The  solution  of  (7)  was  obtained  numerically.  Following  the  scheme 
developed  in  reference  [3],  the  time  interval  t/  was  divided  into  sub¬ 
intervals  At  »  tf/q  and  each  intermediate  time  was  denoted  by  t,  = 
jAt.  Denoting  8j  =  6(lj)  =  0O  +  Rt;,  a,  =  «r(t > )  and  computing  the 
rt -values  of 


ds 

«r(0(s)) 


and  n -values  of  Doy  =  D0(8,)  =  [0.2371  -  0.000574(0;  -  297)]  Hr6  the 
following  expression  was  obtained  for  ot 


D0j  +  1/2 <*;  -  ,)' 

1  _  ..  >=} 


- - Rtj  +  -  Di({;  -  f;— i)'<r;— l 

1  -  v  2 


-  l-  DM  -  2)  L  [(€; -  (m)q  +  <f,  - 

2.  m—  1 


(8) 


In  (8),  W0  -  2)  is  a  unit  step  function  which  vanishes  for  j  <  2.  In 
addition,  it  is  necessary  to  employ  the  initial  conditions  <r()  =  0, 

£o  =  0. 

Calculations  were  performed  on  an  Amdahl  470  computer.  The 
residual  stress  at  the  termination  of  cool-down  is  shown  in  Fig.  1  as 
a  function  of  the  time  t/  allowed  for  the  cooling.  The  reduction  in  the 
residual  stress  with  the  increase  in  cooling-time  t/  is  due  to  viscoelastic 
relaxation.  The  value  of  a  *  39430  kPa  at  t/  =  0  corresponds  to  the 
linearly  elastic  case. 

Additionally  relaxation  of  stresses  occurs,  at  0  *  0/.  at  times  t  >  tf. 
The  calculation  of  these  stresses  requires  a  slight  modification  of 
expressions  (7)  and  (8),  which  is  omitted  here.  The  additonal  relax¬ 
ation,  which  is  drawn  by  the  dashed  lines  in  Fig.  1,  further  reduces  the 
residual  stresses  by  about  3-4  percent. 


Optimal  Cool-Down  Path  for  the  Thermorheologically 
Simple  Model 

It  has  been  shown  previously  |3)  that,  for  a  prescribed  time  1/  for 
cool-down,  there  exists  an  optimal  cool-down  path  0  «  0(f)  which 
minimizes  the  residual  stress  <r/.  A  remarkable  feature  of  this  optimal 
path  iB  that  it  contains  discontinuities  at  times  (  «  0  and  f  *  (/.  It  is 
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where 


J»t/  f  PU  ds  9  _l  . 

1  +  o  I  — —  0(1 'Ml ' 

o  arW 


(13) 


Fig.  2  Relaxation  modu<M£(OoMa<Md&r  approximate  lnvaritona<0(/> 

as  gtvwi  in  squatton  (11)  (daahsd  line)  end  by  exponsnHal  — rt—  determined 
through  coNoceUon  (eoM  line) 


interesting  to  note  that  the  discontinuous  character  of  the  solution 
is  apparently  quite  common,  since  a  similar  result  was  obtained  re¬ 
cently  for  optimal  strain  paths  (5). 

In  order  to  reduce  the  computational  effort  consider  the  approxi¬ 
mate  form  of  the  creep  compliance 

0(1)  =  6(1  +  at9) 

where 

6  =  0.3126  X  10“«  (kPa)~l,  a  =  0.032,  q  =  0.19  (9) 


and 


Also,  in  (13) 


Cl 


3.2  X  lO6® 

I  - 

1  -  v 


<xt  -  exp 


[6480 

0(s) 


-  21.82 


To  determine  the  optimal  path  0(1)  subdivide  t/  into  n  equal  portions 
as  before  and  consider  a  “stair-case”  drop  in  8 


whereby 


and 


0(1)  =  A*,  tk<t<  t*+i‘,  0  <  k  <  n  -  l 


0(1)  =  (A*  -A*-j)  6(1-1*),  1  <k<n 


(14) 


0(1)  =  (Ao-0o)6(t),  6  =  0 
Substitution  of  (14)  into  (13)  gives 

s(t/)  -  A°~9° 


1  + 


a  At  1/orM/)" 
y-o 


n— 1 

+  E 

*-i 


A*  -  A*_ 


■  n  — 1  9 

1  +  a  {At  Z  llar(Aj) 


-(An-,-0/)  (15) 


Equations  (12)  and  (15)  express  the  residual  stress  aft/ )  in  terms 
of  the  unknown  intermediate  temperatures,  Aj  ~  Aft/)  j  =  0, 1, . . . 
n  —  1. 

For  a  stationary  value  it  is  necessary  that  dSft/)/dAj  =  0, ;  =  0, 1 
...  n  —  1.  This  results  in  n  simultaneous  equations  in  the  n  unknown 
Aj  as  follows: 

A0  -  0O  dZ0  A  A,  -  A,-i  dZi  1  1 

k  =  0, 1, ...  n  -  1  (16) 

where,  in  (16)  Z„  -  1,  H(i  -  1)  is  the  unit  step-function  which  van¬ 
ishes  for  negative  arguments,  and 


Equation  (9)  was  derived  from  (5),  with  0  =  75°F  (297°K)  and  when 
Dt,(0)  is  approximated  by  a  constant  value  Do  =  6. 

Furthermore,  it  is  more  convenient  to  employ  the  relaxation 
modulus  Eft)  instead  of  the  compliance  Dft).  This  task  is  accom¬ 
plished  with  the  aid  of  Laplace  transforms.  The  Laplace  transform 
of  Dft )  yields 


pZJ(p)  *■  6 


1  +a 


P9 


pEfp) 


where  Efp)  is  the  Laplace  transform  of  Eft )  (6j.  An  approximate  i 
version,  based  upon  the  expression  (7) 


I  n-1 

Zj  =  1  +  a  At  £  l/«r(A; 

I  y-i 

dZ,  [  n-I 

—*-=  -6480 A*-2 (At  Z  l/«l 
AA*  ( 

The  solution  of  (15)  and  (16)  was  obtained  by  an  iteration  technique.2 
The  resulting  optimal  path  is  shown  in  Fig.  3  for  tf  -  100  min.  It  is 
worth  noting  that  the  value  of  aft/)  is  insensitive  to  the  details  of  the 
“tail-end”  of  the  path,  i.e.,  near  tf.  Cooling  along  the  optimal  path 
results  in  a  residual  stress  aft/)  =  32900  kPa.  This  represents  an  im¬ 
provement  of  about  6  percent  over  the  linear  cool-down  path. 


Eft)* 


1 


pDfp)\p-i/2t 


yields 


Eft)- 


3.2  X  106 
1  +  0.03361°  19 


kPa 


(11) 


The  validity  of  (11)  was  verified  by  comparison  with  a  collocation 
scheme  for  the  inversion  of  ( 10).  The  collocation  method  will  be  dis¬ 
cussed  in  the  following  section  on  nonlinear  effects.  The  comparison 
is  shown  in  Fig.  2. 

Employment  of  (11)  in  (3)  yields 

<r(f/)  -  -C| Sft,)  (12) 


Stress  Effects  on  Relaxation  During  Cool-Down 

We  note  that  the  residual  stresses  reach  levels  which  exceed  half 
the  ultimate  stress  of  epoxy,  au  =  52,200  kPa.  Stresses  of  such  mag¬ 
nitude  were  found  to  influence  the  creep-response  of  epoxy-based 
composites  [9],  The  stress  effect  can  be  incorporated  into  the  "shift 
factor"  «r  of  equation  (4)  as  follows: 


2  We  refrained  from  using  standard  optimization  methods,  because  those 
techniques  presuppose  a  smooth  functional  which  can  be  approximated  by 
quadratic  function  (see.  e.g  .  reference  |8|)  while  the  present  case  revolves  about 
a  discontinuous  function. 
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modal 


Fig.  4  Variation  oi  residual  strata  a  with  lima  I.  Strata  aftacted  relaxation 
(aquation  (17)).  Unaar  (conatant-rata)  cool-down  from  0,  =  438°K  to  8,- 
293°K.  Solid  lina  corraaponda  to  cooling  thna  t,  ~  100  min  and  alto  shows 
ratasatlon  bayond  tarmination  o *  cool-down.  Daahad  Una  corraaponda  to 
t,  =  300  min. 


ar(<r,8) 

\A/8r  —  (B/9r)H[o  —  ac(8)\[o  —  erc(0)]/au{8)  A| 

=  exp|  __  (17) 

In  (17)  8k  is  a  reference  absolute  temperature,  ac  is  a  “threshold 
stress"  beyond  which  stress  starts  to  affect  aT,  <ru  is  the  ultimate 
stress  of  the  resin,  H(  )  is  a  unit  step-function  which  vanishes  for 
negative  arguments  and  A,  B  are  material  constants. 

From  the  data  given  in  reference  [9J,  which  were  obtained  at  a  ref¬ 
erence  temperature  0r  =  293°K,  and  from  recent  data  which  detailed 
the  dependence  of  <ru  on  8  (10),  the  following  values  were  selected: 

A  =  6500,  B  =  2000 

au(6)  •  55200  -  190(0  -  6r) 

<rr(0)  =  13800  -  47.5(0  -  8r). 

where  <ru  and  ac  are  in  kPa. 

A  significant  reduction  in  the  computational  effort  is  achieved  if 
the  relaxation  function  E(t ),  as  given  in  (11),  is  expanded  in  a  sum 
of  exponentials.  To  this  purpose,  the  method  of  collocation  (7]  was 
employed  and  E(t)  was  expressed  as 

£(«)  »  L  £/e'x"  A/  =  10“ 18,  10" 17, . . .  103  (18) 

/-i 

The  values  of  £/,  which  are  not  listed  here,  were  obtained  from 
solving  the  simultaneous  equations 

Z  d  +  X//M-l£/-pE(p)|P-v 

/-i 

with  pE{p)  given  in  (10).  The  form  (18)  of  £(()  is  shown  in  Fig.  2. 
Substitution  of  (17)  and  (18)  in  (3)  and  (4)  yields 

22  /it  d( 

<r(t)  -££/  I  exp|-A,*(t)  +  A, {«')]—  dt'  (19) 
/-i  Jo  at' 

where 


U 


ds 

a7-|0(s),or(s)] 


with  the  form  of  ar  given  in  (17). 

Consider  a  linear  cool-down  from  0O  to  0/  during  a  time  interval  t/. 
Denote  the  constatnt  cooling  rate  by  R  =  (0O  -  0/ ) /f /  and  let 
t'2  =  «£/(!  -  i>).  In  this  case  (19)  reduces  to 


a(t)  -  CiSU) 


where 

£(()  =  ££/  r  exp  [-A/{(t)  +  A/{(t')l<ft'  (20) 

;  Jo 

Divide  the  time  interval  (/  into  n  equal  intervals  At  and  denote  f* 
=  k  At  as  before.  The  initial  condition  is  oo  =  0  and  the  intermediate 
temperatures  are  0*  =  0 o  ~  fit*.  If  we  employ  the  approximation 
that 


ar(<r,0)  «  or(<r*,0*)  for  f*  <  t  <  t*+,  (21) 

then,  as  a  consequence  of  expansion  (18),  it  is  possible  to  integrate  Sit) 
analytically  over  each  subinterval  t*  <  t  <  f*+i-  This  leads  to  a  sub¬ 
stantial  reduction  in  computations. 

Denote  14  =  ar(<r*,7’*)  and  IV*  =  e ~ Vk.  From  (21)  it  follows 
that 

«*-«(!*)=  £  (f,  -  t.-,)IV,-, 

i-l 

Denote  further 

Z  Urn  -  tm-|)W4,-|  -  £*  (r=  1,  ..k.  l,..n) 

mm0 

and 


<t>rkm*rk+  Wr-,(tr  -tr— ,) 
Straightforward  manipulations  lead  to 


1 
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Fig.  S  Two-st«f>  cool-down;  variation  of  tho  rooktoa)  strns  a  wHh  tho  Ui- 
tormodtolo  lavol  twnpcrituri  00  =  4W°K,  6,  =  M3°K,  f,  =  100  min.  Not* 
minimal  realduat  *f teas  at  0,  =  333 °K. 


22  (  k 

S*  -S(t*)  =  E  El  E  (A/ltV-,)-' 

/-I  lr*l 

X  [exp  (X/<^>rfc)  -  exp  (A/^rktij  (22) 

and 

Ok  =  c2S* 

For  times  beyond  termination  of  cool-down  f*  >  tf,  the  temperature 
retains  a  constant  value  0/  but  the  stress  <r«  continues  to  relax.  The 
computational  scheme  (22)  remains  valid  except  that  now,  «r  (a*.®* ) 
«  ar(tr*,0/),  (tK  >  tf). 

Results  for  a  constant  cooling  rate  from  do  =  438°K  down  to  6/  = 
293°K  are  shown  in  Fig.  4  for  tf  ~  100  min,  followed  by  a  relaxation 
at  293°  K,  and  for  t /  =  300  min.  Note  that  a  prolongation  in  the  cool¬ 
down  time  does  not  effectively  reduce  the  residual  stresses.  Com¬ 
parison  with  results  based  upon  the  lineal  thermorheologically  com¬ 
plex  model  shows  that  when  stress  effects  are  incorporated  into  «r. 
equation  (17),  the  relaxation  is  enhanced  and  the  residual  stresses 
reduce  by  an  additional  amount  of  10  percent. 

Finally,  consider  the  case  of  a  “two-step  cooling."  Tn  this  case  the 
temperature  is  dropped  suddenly  from  0O  to  0;„  held  at  0  =  0t,  for  a 
period  t  =  </  and  dropped  sharply  again  to  0  =  0/  at  t  =  f/.  For  such 
cooling-history  equation  (19)  yields 

=  cmS(tf)  +  <yE  Ef  (23) 

l 

In  equation  (23) 

Cm  =  -  (00  —  9b),  Cf  =  - -  (0),  -  6f ) 

1  -  V  1  -  l> 

and 

Si  =  E 
l 

where  &  is  defined  and  evaluated  as  before. 

Following  a  numerical  scheme  similar  to  (22),  the  dependence  of 
the  residual  stress  «r(f/)  on  the  intermediate  temperature  0*  has  been 


computed  for  a  “hold-time”  tf m  100  min.  In  the  calculations  we  took 
0o  =  438®  K  and  0/  -  293°  K.  The  results  are  shown  in  Fig.  5. 

Note  that  the  residual  stress  attains  a  minimum  at  the  intermediate 
temperature  0t  -  333°K.  A  similar  behavior  was  observed  experi¬ 
mentally  for  epoxy  coupons  reinforced  by  unidirectional  fibers.  Those 
coupons  were  cooled -down  in  two-steps  and  a  qualitative  assessment 
of  the  residual  stresses  was  obtained  by  means  of  a  frozen  fringe 
pattern.  It  was  observed  that  a  minimal  fringe  count  was  attained  at 
an  intermediate  value  of  temperature  0(,(0/  >  0s  >  0o)  3 

Concluding  Remarks 

It  has  been  shown  that  a  linear  elastic  analysis  of  residual  stresses 
overestimates  their  magnitude  by  more  than  20  percent.  If  account 
is  taken  of  the  viscoelastic  response  of  the  material  it  may  be  possible 
to  design  a  cool-down  path  that  will  yield  a  reduced  value  of  residual 
stresses.  It  seems  that  the  details  of  the  cool-down  path  influence  the 
magnitude  of  the  residual  stresses  more  noticeably  than  the  duration 
of  the  cooling  period. 

The  computations  conducted  in  the  present  paper  are  based  upon 
incomplete  data.  It  is  intended  to  modify  the  calculations  and  analysis 
in  light  of  new  data  as  they  become  available.  It  should  be  noted  that 
the  effects  of  aging  were  not  considered  in  the  present  formulation. 
These  effects  may  be  of  prime  importance  because  it  is  conceivable 
that  the  newly  created  epoxy  may  lie  undergoing  rapid  aging  while 
it  is  being  cooled  down.  An  inquiry  into  this  matter  is  now  under 
way. 
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Extended  Abstract 

Increasing  concern  with  the  effects  of  temperature  and  moisture 
on  the  performance  of  composites  demands  the  development  of  accelerated 
environmental  conditioning  techniques.  These  conditioning  schemes 
attempt  to  simulate  extreme  exposures  to  moisture  and  temperature  which 
may  be  incurred  during  service  conditions  of  the  composite  material. 

Since  the  diffusion  of  moisture  proceeds  at  extremely  slow  rates, 
and  may  require  months  or  years  until  saturation  is  reached  at  room 
temperature,  the  conditioning  process  is  often  accelerated  by  the 
employment  of  elevated  conditioning  temperatures  which  speed  up  the 
rate  of  moisture  sorption  and  shorten  the  time  of  the  experimental 
task. 

It  was  observed  that  certain  conditioning  schemes  may  result  in 
microcracking  and  damage  to  the  composite  laminates.  In  order  to 
develop  a  time-saving,  yet  non-damaging,  environmental  conditioning 
method,  a  need  exists  for  evaluating  the  internal  stresses  which  arise 
due  to  moisture  and  temperature  and  assess  the  relative  merits  of 
various  conditioning  schemes. 

In  this  paper,  the  moisture  and  temperature  effects  are  formulated 
by  means  of  a  viscoelastic  response.  Elastic  solutions  are  provided  for 
purpose  of  comparison.  The  computations  incorporate  realistic  material 
parameters  that  are  based  on  recent  data  and  results  are  presented  for 
symmetric  laminates  containing  0°  and  90°  plies.  The  initial  stresses 
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due  to  cool-down  from  cure  and  due  to  moisture  sorption  during  storage 
prior  to  conditioning  are  calculated  by  means  of  elastic  analysis. 

Since  heat  diffuses  extremely  rapidly  and  equilibrates  throughout 
the  laminate  during  times  that  are  much  shorter  than  all  other 
characteristic  times,  we  discard  the  process  of  heat  conduction  and 
consider  spatially  uniform  temperature  profiles  for  the  entire  condition¬ 
ing  stage.  On  the  other  hand,  we  account  for  the  predominant  effect  of 
temperature  on  mositure  diffusivity  in  solving  for  the  moisture  diffusion 
which,  due  to  its  extremely  slow  progress,  is  described  by  a  time- 
transient  and  spatially  non-uniform  process  throughout  the  conditioning 
stage. 

Moisture  and  temperature  effect  the  response  of  composites  in  two 
counter-acting  ways.  On  one  hand,  both  factors  induce  swelling  strains 
into  the  material,  particularly  into  the  matrix  component  of  the 
composite.  Since  the  natural  tendency  to  swell  is  resisted  by 
geometric  constraints  which  are  always  present  in  structures— and 
especially  in  cross-ply  laminates— both  moisture  and  temperature  can 
be  viewed  as  direct  causes  for  internal  stresses  and  are  recognized 
as  such  in  both  the  elastic  and  viscoelastic  analyses. 

On  the  ohter  hand,  both  moisture  and  temperature  strongly  affect 
some  of  the  material  moduli  of  fiber-reinforced  composites.  These 
effects  are  most  pronounced  in  the  softening  and  accelerated  rates  of 
creep  (or  relaxation)  that  occur  with  rising  temperatures  and  moisture 
levels.  These  last-mentioned  effects  tend  to  counteract  the  afore¬ 
mentioned  influences,  so  that  with  increase  in  temperature  and  moisture 
an  actual  competition  takes  place  between  stressing  on  one  hand  and 
stress  relief  on  the  other.  While  this  "competition"  is  accounted 
for  in  the  viscoelastic  analysis,  it  cannot  be  considered  within  the 
context  of  the  elastic  response. 

All  the  above  mentioned  factors  are  included  and  featured  in  our 
viscoelastic  analysis.  In  all  our  computations,  we  disregard  edge 
effects  and  consider  only  those  stress  and  strain  fields  which  are 
several  ply-thickness  distance  away  from  the  edges.  The  elastic  and 
viscoelastic  solutions  are  evaluated  up  to  the  times  at  which  the 


laminates  reach  moisture  saturation,  whereupon  they  are  cooled  back 
to  room  temperature. 

Some  representative  results  are  shown  in  Figures  1  and  2.  Both 
figures  show  the  stress  a  which  arises  in  an  8-ply  laminate,  which 
consists  of  two  outer  0°  plies  and  two  inner  90°  plies  arranged 
symmetrically  about  the  middle  plane.  The  y  axis  coincides  with  the 
longitudinal  (parallel  to  fibers)  direction  of  the  90°  plies  and 
corresponds  to  the  transverse  direction  in  the  outer,  0°  plies.  The 
stress  profiles  are  plotted  vs.  the  depth  coordinate  z  across  the 
thickness  of  the  laminate  (z  =  0  at  the  middle  plane)  at  times 
t  =  1,  100,  1000  and  10,000  minutes  into  the  conditioning  stage  and 
also  at  saturation,  just  prior  and  just  past  cooling  back  to  room 
temperature. 

Figure  1  provides  a  comparison  between  two  conditioning  schemes. 
Exposure  to  98%  R.H.  at  180°  F  (solid  line),  for  which  saturation  is 
reached  after  22  days,  and  exposure  to  98%  R.H.  at  150°  F  (dashed  line), 
for  which  saturation  is  attained  after  45  days.  Note  the  substantial 
differences  in  the  stress  profiles  for  the  two  conditioning  schemes. 

In  particular,  upon  return  to  room  temperature  the  stresses  described 
by  the  solid  lines  are  of  opposite  sign,  and  somewhat  larger,  than 
those  given  by  dashed  lines.  In  all  other  lay-ups,  for  which  computa¬ 
tions  were  performed,  it  was  similarly  found  that  conditioning  at 
150°  F  resulted  in  lower  values  for  the  final  stresses. 

Figure  2  provides  a  comparison  between  viscoelastic  values 
(shown  in  solid  lines)  and  elastic  results  (dashed  lines).  The 
computations  were  performed  for  conditioning  at  98%  R.H.  and  180°  F. 

Note  the  vast  discrepancy  between  the  two  sets  of  results,  especially 
with  increasing  time,  indicating  that  by  discarding  stress  relaxation 
the  elastic  analysis  yields  results  that  are  substantially  off  the 
mark. 


Thermal  Buckling  of  Free  Plates 
due  to 

Non-Uniform  Temperature  Distribution 


by 

Y.  Weitsman* 


This  investigation  concerns  the  buckling  of  thin  isotropic  plates  due 
to  uneven  temperature  distribution  over  their  surface. 

The  problem  is  motivated  by  the  processing  and  operational  conditions 
of  fiber  reinforced  composite  laminates.  It  was  noticed  that  during  the 
cool-down  of  large  composite  structures  temperature  variations  of  up  to 
40°F  could  occur  in  several  spotty  locations.  In  addition,  localized  "hot 
spots"  may  be  anticipated  in  various  operational  circumstances. 

In  general,  the  actual  determination  of  buckling  conditions  requires 
the  employment  of  approximate  methods  of  analysis  and  computation.  However, 
the  success  of  any  approximate  technique  depends  crucially  on  the  selection 
of  reasonable  "shape  functions"  which  reflect  the  buckled  shape  of  the 
structure. 

In  this  Note  we  consider  an  idealized,  rotationally  symmetric  case 
which  is  amenable  to  an  exact  analytical  solution.  This  solution  offers  a 
model  for  a  suitable  "shape  function”  and  also  provides  an  assessment  for 
the  circumstances  under  which  thermal  buckling  may  become  a  significant  de¬ 
sign  factor. 

Analysis 

Consider  a  thin  isotropic,  elastic  plate  of  thickness  h  and  radius  b, 
which  is  traction-free  at  r  *  b. 

Let  E,  v  and  a  denote  the  Young's  Modulus,  Poisson's  ratio  and  the 
coefficient  of  thermal  expansion,  respectively.  The  bending  stiffness  of 
the  plate  is  given  by  D  =  Eh^/1 2(1 -v^) .  Let  r  denote  the  radial  coordinate 
and  consider  the  following  temperature  distribution 

T(r)  -  yi-rV|  *  Tb  (1) 

^Professor,  Mechanics  &  Materials  Research  Center,  Civil  Engineering 
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Elementary  computations  [1]  yield  the  following  expression  for  the  radial 
stress  op: 

0r  -  l  Ec  T0[-l+(r/b)2]  (2) 

Note  that  or  <  0  throughout  the  plate. 

Let  <f>  designate  the  slope  of  the  buckled  plate.  Employing  (2),  we  ob¬ 
tain  the  following  expression  for  41  [2]: 

rV  +  r<f>'  -  [{^)2r4-k2r2  +  1]  *  =  0  (3) 


where  k 
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EoT  h 
0 

4D 


=  3(l-v2)  — 
fr 


,  and  primes  designate  derivatives 


with  respect  to  r.  It  is  interesting  to  note  that  (3)  does  not  depend  on 

Young's  modulus.  This  is  due  to  the  fact  that  E  contributes  to  both  the 

stress  o  and  the  stiffness  D. 
r 

2  k 

Introduce  the  new  variables  p  =  r  ,  z  =  ^  p  and  n. 

Let  <t>(r)  =  exp  (-  |j  p^  n(z).  Subsitution  into  (3)  gives 

zn"  +  ( 2-z ) n '  -  (1  -  ^)n  =  0  (4) 


where,  in  (4),  primes  indicate  derivatives  with  respect  to  z. 


Equation  (4)  is  Kummer's  equation  [3]  [4]  which  has,  in  general,  two 

independent  solutions.  In  the  present  case,  one  solution  contains  the 

term  log  z.  This  logarithmic  solution  is  rejected  because  its  leading 

-1  k  2 

term  leads  to  the  expression  $  =  Ag[r  exp  (-  ^  r  )  + - ]  which  is 

singular  at  r  =  0. 


The  remaining  solution  yields,  after  several  straightforward  manipu¬ 
lations 


♦  *  Ai  x  exp  (-  x2)  U(1  -  ,  2,  kbx2) 


In  (5),  U(a,  y ,  z) 


=  rfv) r(a+n)  zn 
r(a)  r(Y+n)  rTT 


n=o 


x  =  r/b  and  Aj  is  an  arbitrary  constant. 


(5) 


For  a  traction-free  plate  both  bending  moment  Mf  and  shear  force 
Qr  vanish  at  r  =  b.  However,  since  Qr  =  Nr<*>  =  ho^  it  follows  from 
(2)  that  Qr(b)  =  0.  The  remaining  condition,  Mr(b)  =  0,  gives 


Substitution  of  (5)  into  (6)  yields,  after  several  intermediate  cal¬ 
culations,  the  following  transcendental  equation  in  x 

(1+v-A)  U(1  -  X/4 ,  2, A)  +  x(l-A/4)  U ( 2  -  X/4,  3,X)  =  0  (7) 

where,  in  (7) 

x  =  kb  =  V30-v2)aT  (b/h)  (8) 

w«  0 

The  smallest  non-zero  root  of  (7)  determines  the  critical  temperature  Tq 

which  would  cause  thermal  buckling.  Substituting  this  value  of  A  =  kb 

dU 

into  (5)  we  get  the  slope  of  the  buckled  shaped  -  ^  .  The  deflected 
shape  W  is  obtained  by  numerical  integration  of  (5). 

Computations  were  performed  with  typical  values  for  quasi-isotropic 

-6  ° 

graphite/epoxy  and  glass/epoxy  laminates,  namely  a  =  7  x  10"  in/in/1  F, 
v  =  0.3,  and  with  T  =  20°F.  The  smallest  root  of  equation  (7)  was  found 

0  L 

to  be  xCR  =  2.405.  Substitution  in  (8)  yields  ^  =  124.3 

It  follows  that  unstiffened  quasi-isotropic  laminates  with  in-plane 
dimensions  that  exceed  100  times  the  thickness  may  be  susceptible  to  thermal 
buckling  if  the  non-uniformity  in  the  temperature  distribution  exceeds 
20°F. 

A  typical  buckled  shape  is  shown  in  Fig.  1,  where  the  non-dimension- 
alized  deflection  w  =  W/b  is  plotted  vs.  the  non-dimensional  radial  distance 
x  a  r/b.  As  in  all  buckling  problems,  the  amplitude  of  the  deflection 
remains  indeterminate. 
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Figure  Title 

Shape  of  Deflection  w(x)  vs.  Dimensionless  Distance  x  =  £withv=  0.3 
The  Critical  Value  of  A  is  Aj.r  =  2.405. 
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A  METHOD  FOR  DETERMINING  THE  MODE  I  DELAMINATION  FRACTURE  TOUGHNESS 
OF  ELASTIC  AND  VISCOELASTIC  COMPOSITE  MATERIALS* 


by 

D.  F.  Devitt**,  R.  A.  Schapery+  and  W.  L.  Bradley++ 

ABSTRACT 

A  method  using  a  double  cantilever,  laminated  beam  specimen 
is  developed  to  obtain  delamination  fracture  toughness  as  a  function 
of  the  rate  of  crack  propagation.  First,  the  relation  between  energy 
release  rate,  applied  displacement,  and  various  laminate  parameters  is 
derived  using  a  large  displacement,  small  strain  theory.  Experiments 
employing  unidirectional  glass/epoxy  composites  of  three  different  thick¬ 
nesses  and  a  wide  range  of  loading  rates  are  then  described.  Although  the 
beam  deflections  and  rotations  are  very  large,  good  agreement  between 

measured  and  predicted  beam  compliance  is  demonstrated.  The  energy  release 

•  *0  1 

rate  Gc  and  crack  speed  L  are  shown  to  obey  the  power  law  Gc  ~  L  ; 

essentially  the  same  result  is  obtained  for  ail  three  laminate  thicknesses. 
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INTRODUCTION 


Predictions  of  the  useful  life  for  composite  material  components 
require  a  comprehensive  understanding  of  the  materials  response  to 
complex  load  histories  in  anticipated  service  environments  in  the 
presence  of  defects  such  as  microcracks,  voids,  and  delaminations.  Such 
defects  can  occur  during  manufacturing  or  may  develop  in  service, 
causing  structural  degradation  or  failures  at  stresses  well  below  the 
strength  levels  expected  for  defect  free  material.  Linear  elastic 
fracture  mechanics  (LEFM)  has  been  developed  to  deal  with  crack-like 
defects  by  relating  defect  geometry  and  design  stress  to  a  material 
response,  normally  called  the  fracture  toughness  of  the  material.  The 
fracture  toughness  of  a  material  is  usually  characterized  by  critical 
energy  release  rate  (G  )  or  the  critical  stress  intensity  K  . 

In  an  ideal  monolithic  isotropic  material  the  fracture  toughness  is 
independent  of  the  orientation  of  the  crack  plane  for  a  given  mode  of 
deformation  at  the  crack  tip,  such  as  the  opening  mode.  However,  in 
an  anisotropic  composite  material,  the  material  response  may  vary 
considerably  depending  on  the  plane  of  fracture  and  the  energy 
dissipative  processes  involved.  Work  to  date  has  indicated  toughness 
values  of  14-28  ksi/m  for  center  notched  tensile  specimens  of  fiber- 
reinforced  plastic  laminates  [1],  which  compare  favorably  with  aluminum 
alloys  often  used  in  the  aircraft  industry  which  have  Klc  values  of 
21-40  ksi/Tn"  [2]. 

Growth  of  interlaminar  flaws  (delamination)  is  an  important  part 


of  the  failure  process  in  many  laminates  [3,  4].  Compressive  fatigue 
appears  to  be  an  especially  severe  type  of  loading  in  producing 
delaminations  [5,  6];  out  of  plane  stresses  developed  through  compressive 
loading  and  local  buckling  are  thought  to  be  the  primary  cause  of  such 
delamination  type  fractures.  These  observations  indicate  that  the 
delamination  fracture  toughness  may  be  the  critical  toughness  parameter 
for  fatigue  stressing  where  the  minimum  stresses  are  compressive. 

While  considerable  effort  has  been  expended  to  define  fracture 
toughness  for  tensile  loading  of  laminates  with  flaws  normal  to  lamina 
planes,  very  little  has  been  done  to  better  define  and  understand  the 
delamination  fracture  behavior.  Apparently,  the  few  studies  that  have 
been  conducted  to  characterize  delamination  have  utilized  specimen 
geometries  that  give  a  mixed  tensile  and  shear  stress  state  at  the 
delamination  crack  tip  [3,  7,  8]  (see  Figure  la):  precluded,  therefore, 
is  determination  of  Klc  or  G-|c  for  the  opening  mode  of  delamination  in 
which  only  tensile  stresses  exist  at  the  crack  tip.  In  order  to  fully 
characterize  delamination  fracture  toughness  we  believe  it  is  necessary 
to  study  the  effect  of  various  proportions  of  tensile  and  shearing 
stresses  at  the  crack  tip,  including  the  opening  mode. 

The  objective  of  the  investigation  described  herein  has  been  to 
develop  an  experimental  approach  with  the  associated  analysis  to  obtain 
the  critical  energy  release  rate  for  the  opening  mode  of  delamination. 

A  double  cantilever  beam  geometry  was  chosen  to  give  essentially 
pure  opening  mode  fracture  (Figure  lb).  To  allow  the  study  of  thin 
laminates,  nonlinear  beam  theory  was  used  in  the  analysis.  Experimental 
measurements  were  made  on  Scotchply  (a  composite  consisting  of  epoxy 
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matrix  reinforced  with  continuous  glass  fibers)  to  verify  the  analysis 
and  allow  determination  of  the  delamination  fracture  toughness  of 
this  material.  The  analytical  approach  will  be  presented  in  the  next 
section,  followed  by  a  description  of  the  experimental  program.  The 
results,  including  raw  data  and  reduced  data  will  then  be  presented  with 
appropriated  discussion,  from  which  final  conclusions  will  be  drawn. 

The  thin,  double  cantilever  beam  geometry  gives  rise  to  stable 
crack  growth,  and  therefore  is  particularly  suited  to  determine  the 
relation  between  slow  crack  speed  and  the  fracture  toughness.  A 
portion  of  our  study  therefore  has  been  devoted  to  characterizing 
delamination  growth  rate  in  terms  of  energy  release  rate. 

ANALYTICAL  PROCEDURES 

Consider  the  beam  specimen  in  Figure  lb,  in  which  the  delamination 
crack  tip  is  at  the  point  x  =  a.  The  energy  release  rate  associated 
with  a  virtual  crack  growth  of  an  amount  da  is,  by  definition,  the 
mechanical  energy  that  becomes  available  at  the  crack  tip  per  unit  area 
of  new  surface.  In  terms  of  the  total  strain  energy  in  an  elastic  beam, 
W,  this  release  rate  is 

<’> 

where  the  derivative  is  evaluated  for  constant  beam  tip  displacement  a. 
Also,  B  is  specimen  width  normal  to  the  page.  The  value  of  G  at  which 
the  crack  actually  starts  to  propagate  is  the  critical  energy  release 
rate  or  fracture  toughness,  Gc> 

Evaluation  of  Equation  (1)  for  our  specimen  geometry  (Figure  2) 


*  ft 


requires  deteriiiination  of  the  strain  energy  stored  in  the  double 
cantilevered  beam  as  a  function  of  the  parameters  of  grip  displace¬ 
ment  (a),  moment  arm  L-d,  moment  of  inertia  I,  and  the  modulus  of 
elasticity,  E,  in  the  length  direction  of  the  beam.  The  calculation 
of  the  strain  energy  stored  in  a  nonlinear  beam  in  terms  of  these 
parameters  will  be  presented  next. 


Strain  Energy  in  a  Non! inear  Beam 

The  strain  energy  in  a  cantilever  beam  under  a  concentrated  end 

load  is  found  from  elementary  beam  theory  to  be, 

L  „ 
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where  M  is  the  bending  moment,  E  is  the  axial  modulus  of  elasticity,  I  is 
the  area  moment  of  inertia,  and  L  is  the  underformed  beam  length. 

As  a  result  of  the  low  flexural  rigidity  of  thin  composite  specimens, 
large  deflections  and  rotations  may  be  present  during  the  test  and  need 
to  be  accounted  for  in  the  analysis.  Such  a  correction  to  linear  beam 
theory  has  been  made  by  Bisshopp  and  Drucker  [10],  where  they  included 
the  square  of  the  rotation  in  the  strain-displacement  relation  and  tiie 
effect  of  the  beam  shortening  of  the  moment  arm  due  to  large  deflections. 

A  linear  stress-strain  equation  and  small  strains  were  assumed. 

The  nomenclature  and  geometric  relationships  used  by  Bisshopp  and 
Drucker  are  shown  in  Figure  2.  The  principal  results  of  their  analysis 
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(^)?  QI  0  {sin4>0-  sin4>)  d<j> 


(3) 


_  /EK?  ,  0 


A  =  (tJo-) 


W 


sjn$ 


{sin<}>0  -  sin<t>) 


17? 


d<j> 


(4) 


where  4>0  is  the  angle  of  slope  at  the  loaded  end  of  the  beam,  <j>  is  the 
slope  of  the  beam  at  intermediate  points,  P  is  the  load,  and  the  remaining 
terms  were  defined  previously.  After  determining  the  appropriate  change  of 
variables  and  transformations  Bisshopp  and  Drucker  were  able  to  rearrange 
Equations  (3)  and  (4)  into  canonical  elliptic  form  as  shown  below: 

/PL2x  1/2  =  F(k)  -  F(k,  *)  (5) 

'El  1 
and 


f  -  1  -  2  ( 


E< 

X 

1  -  ♦) 

~n 

1  -  f(k,  <t> 

(6) 


where  F(k)  and  E(k)  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind  respectively,  F(k,  <f> )  and  E(k,  <j>)  are  incomplete  elliptic 
integrals  of  the  second  kind,  and  k  and  $  are  elliptic  parameters  [11]. 

A  PL2  V2 

Since  the  results  for  j-  and  (^-)  are  single-valued  functions  of  the 
same  elliptic  parameters,  they  are  implicitly  related  through  Equations  (5) 
and  (6),  with  relationship  presented  in  Table  1  and  Figure  3. 

Next,  the  strain  energy  stored  in  the  two  cantilevered  beams  must 
be  calculated.  Combining  the  general  definition  for  strain  energy  W-j  in 
a  beam  (Equation  2)  with  the  specific  result  from  non-linear  beam  analysis 
(Equation  3),  one  may  express  the  total  strain  energy  W  *  2W^  in  terms  of 
elliptical  integrals  as  follows: 


i 


6 


KT  ■  <F(k)  -  F(k,  *))2  (  ffij-Hfc-*]  *  2(k2  -  1)  )  ID 

The  results  of  such  a  calculation  are  presented  in  Figure  4  and 
Table  1.  Next,  the  energy  release  rate  as  defined  in  Equation  1  will 
be  evaluated  using  the  strain  energy  for  a  nonlinear  beam  as  calculated 
in  Equation  8. 


Energy  Release  Rate  (G)  for  a  Nonlinear  Beam 


It  is  clear  from  Figure  4  that  increasing  the  load  (or  grip) 
displacement  "A"  for  a  fixed  crack  length  "L"  will  increase  the  strain 
energy  stored  in  the  specimen,  "W".  From  the  same  Figure  4  it  is  also 
clear  that  crack  extension  (increasing  L)  which  occurs  at  constant  grip 
displacement  (a)  results  in  a  decrease  in  the  strain  energy  (W)  stored 
in  the  specimen,  as  shown  schematically  in  Figure  5.  As  the  specimen 
with  length  (or  crack  size)  L-j  is  loaded  from  a-j  to  A^,  the  strain 
energy  W  is  seen  to  increase.  The  strain  energy  release  rate 
AW 


Ga  -  ■  BaL 


(8) 


is  also  seen  to  increase  as  the  test  goes  from  a-b.  Suppose  at  "b"  the 
critical  energy  release  rate  is  achieved  and  crack  extension  occurs, 
changing  the  crack  length  from  L-j  to  At  point  "c"  the  energy 
release  rate  is  again  insufficient  to  give  crack  extension,  so  that 
further  loading  to  point  "d"  is  necessary,  which  again  corresponds  to 
G  equal  to  the  critical  energy  release  rate.  Crack  extension  occurs  from 
Lg  to  Lj  as  one  moves  along  d-e.  In  an  actual  test  the  loading  and  crack 
extension  occur  more  or  less  continuously,  instead  of  in  steps.  Indeed  the 
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energy  release  rate  may  be  calculated  as  --g-  for  constant  A.  The 
path  a-c-e  for  G  =  Gc  represented  by  the  dotted  line  in  Figure  5  would 
be  followed  in  an  actual  experiment. 

The  present  theoretical  prediction  of  G  can  be  used  for  viscoelastic 
materials  if  the  axial  stress-strain  relation  is  elastic;  this  is  a  good 
approximation  for  our  tests  since  we  have  employed  only  unidirectional 
laminates  with  fibers  parallel  to  the  beam  axis.  However,  the  resin  has 
a  significant  influence  on  the  energy  absorbed  at  the  crack  tip,  and 
therefore  the  beam  will  be  locally  viscoelastic.  Consequently,  crack 
growth  rate  L  depends  on  the  stress  intensity  factor,  Kj  and  therefore 
on  the  critical  energy  release  rate  [12].  For  the  usual  case  in  which 
dL/dG  >  0,  a  slower  crosshead  speed  results  in  a  slower  crack  growth 
rate;  thus,  a  different  path  with  smaller  values  of  G  would  be  traversed 

as  shown  by  the  dotted  line  in  Figure  5. 

WL 

Since  is  a  single-valued  function  of  a/L,  as  seen  in  Figure  4, 
one  may  write 


=  f(A/L) 


Differentiation  of  "W"  at  constant  a  gives 


1  3W 
1  3L 


(9) 

(10) 


Next,  define  S’, 


S  =  d  1 og( WL/E I )/d  log(A/L) 
Therefore,  from  Equations  (9)  -  (11), 

r  -  ^W  -  E 1  /WL  \  / 1  ,  q’x 

G  "-Bar-  {rr}  (1  +  s) 


(id 


(12) 
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The  logarithmic  slope  ?  is  drawn  in  Figure  6. 

A  rearrangement  of  Equation  12  yields  a  nondimensional  form  of  the 
energy  release  rate  for  one  of  the  two  beams. 


1 
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GBL2  _  I  WL 
El  ‘  2  El 


(1  +  S) 


(13) 


_  2 
Since  both  WL/EI  and  S  in  Equation  13  are  functions  of  A/L,  GBL  /El  is  an 

implicit  function  of  A/L.  Using  Equation  13  in  conjunction  with 

2 

Figures  4  and  6,  GBL  /El  as  a  function  of  a/L  has  been  evaluated  and 
presented  graphically  in  Figure  7.  Suppose  that  during  quasi-static 
crack  growth,  one  determines  A/L  at  any  given  time  in  a  delimination 
fracture  test:  Figure  7  then  allows  one  to  determine  Gc  directly  as 
G  =  G  for  a  growing  crack. 

In  summary,  the  relationshops  between  a/L  and  P,  W  and  G  have  been 
determined  in  terms  of  elliptical  integrals  using  nonlinear  beam 
analysis.  These  results  are  summarized  in  Figures  3,  4  and  7  and 
Table  1  and  will  be  used  to  analyze  the  experimental  results  to  be 
described  in  the  next  section. 


EXPERIMENTAL  PROCEDURES  AND  RESULTS 


The  experimental  program  to  be  described  in  this  section  had  two 
objectives:  (1)  to  experimentally  verify  the  nonlinear  beam  analysis; 
and  (2)  to  measure  the  critical  energy  release  rate  as  a  function  of 
crack  growth  rate,  as  suggested  by  Schapery  [12],  In  this  section, 
specimen  preparation  and  testing  will  be  described,  typical  raw  data 
will  be  presented,  and  reduction  of  the  data  to  critical  energy  release 
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rate  values  (Gc)  using  the  previously  described  nonlinear  beam  analysis 
will  be  detailed.  Finally,  a  summary  of  all  of  the  results  will  be 
presented  in  tabular  and  graphical  form. 

Specimen  Preparation 

Scotchply,  which  is  an  E-glass  reinforced  type  1003  epoxy,  was 
selected  for  this  study  because  it  is  translucent,  which  makes  the 
visual  measurement  of  the  moving  crack  front  during  testing  quite  easy. 
Eight,  twelve,  and  sixteen  ply  12"  x  12"  panels  were  laid  up  using 
prepreg  tape  containing  continuous  fiberglass  filaments.  All  lamina 
had  the  same  fiber  orientation  in  the  laminates.  After  curing,  the 
12"  square  laminates  were  cut  parallel  to  the  fiber  direction  into  one 
inch  wide  unidirectional  test  strips.  The  cured  lamina  thickness  is 
about  .00832  inches,  giving  laminate  thicknesses  of  0.067,  0.102  and  0.133 
inches,  respectively,  for  the  8,  12  and  16  ply  specimens. 

Test  Procedure 

The  specimens  were  tested  in  an  ambient  environment  (approximately 
75° F  and  50%  RH)  using  an  Instron  tensile  test  machine  with  a  one 
thousand  pound  load  cell;  a  twenty  pound  full  scale  range  was  used. 

Special  grips  were  designed  to  allow  the  load  to  be  applied  along  a 
fixed  line  of  action  while  permitting  a  virtually  free  rotation  of  the 
loaded  ends  of  the  double  cantilever  beam,  as  was  assumed  in  the 
analysis  (see  Figure  2).  Load  was  monitored  as  a  function  of  time  using 
a  strip  chart  recorder.  The  displacement  versus  time  was  calculated  from 
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the  crosshead  speed.  The  crack  position  as  a  function  of  time  was  noted 
visually  as  the  crack  front  passed  reference  marks  on  the  specimens.  A 
specimen  test  in  progress  is  shown  in  Figure  8.  Typical  experimental 
results  are  presented  in  Table  2  for  a  twelve  ply  specimen  tested  at  a 
crosshead  speed  of  1  inch/minute.  Five  each  of  the  eight,  twelve 
and  sixteen  ply  specimens  were  tested  at  a  crosshead  speed  at  one  inch/ 
minute.  One  specimen  was  tested  with  a  crosshead  speed  that  vanned 
from  .01  inch/minute  to  10  inch/minute. 

Data  Reduction 

The  measured  data  included  load  (P),  grip  displacement  (2a),  and 
crack  length  (1);  they  can  be  easily  reduced  to  give  the  critical  energy 
release  rate  (Gc)  using  the  analysis  described  in  an  earlier  section, 
especially  the  results  of  that  analysis  tabulated  in  Table  1.  After 

I 

determining  A/L,  Lagrangian  interpolation  was  used  to  calculate 

2  — 

PL  /El,  WL/EI,  and  S.  The  results  of  these  interpolations,  along  with  the 
subsequent  calculation  of  the  critical  energy  release  rate  Gc  using 
Equation  13,  are  summarized  in  Table  2.  The  crack  growth  rate  was 
calculated  from  the  measured  values  of  grip  displacement  versus  crack 
extension  (L),  knowing  the  rate  of  grip  displacement.  From  the  measured 
values  of  P,  2a,  and  L,  the  effective  flexural  rigidity  term  was 
determined. 

Experimental  Results 

Results  taken  from  8,  12,  and  16  ply  specimens  have  been  superimposed 

PL2 

on  Figure  4,  as  seen  in  Figure  9.  While  can  be  predicted  directly  by 


the  analysis  as  shown  by  the  solid  line,  the  values  may  also  be  determined 
directly  from  measured  quantities.  The  accuracy  of  the  analysis  is 
clearly  demonstrated  by  the  fact  that  the  experimental  data  points  for 
all  three  specimen  thickness  fall  nicely  on  the  curve  predicted  from  the 
analysis.  Another  way  of  demonstrating  the  accuracy  of  the  analysis  is 
to  predict  the  crack  length  L  at  various  times  in  a  test  using  the 
measured  values  of  P  and  2a.  The  predicted  values  are  compared  with  the 
experimental  values  in  Table  3  for  a  twelve  ply  specimen.  The  agreement 
is  seen  to  be  excellent. 

The  critical  energy  release  rate  results  are  presented  graphically 
in  Figure  10  as  a  function  of  crack  growth  rate,  L.  As  expected,  higher 
crack  growth  rates  required  larger  critical  energy  release  rates.  Further¬ 
more,  the  sixteen  ply  specimens  had  a  higher  Gc  value  than  did  the  eight 
ply  specimen,  though  only  slightly.  The  modulus  of  elasticity  values 
determined  in  the  analysis  indicated  values  of  6.35  x  106,  and 
4.95  x  106  psi  for  the  8,12,  and  16  ply  specimens.  This  variation  seems 
to  reflect  a  decreasing  fiber  volume  fraction  with  laminate  thickness, 
which  was  subsequently  confirmed  using  density  measurements.  This 
type  of  variation  is  reasonable  as  it  is  more  difficult  to  bleed  out  the 
excess  resin  in  the  thicker  laminates,  giving  them  a  higher  volume 
fraction  of  residual  resin,  and  therefore,  a  lower  modulus  of  elasticity. 
The  higher  Gc  values  observed  for  the  thicker  laminates  also  suggest  that 
a  greater  volume  fraction  of  resin  may  slightly  enhance  the  fracture 
toughness.  Photomicrographs  of  the  three  specimens  are  shown  in 
Figures  lla-lld.  The  fiber  distribution  in  all  three  specimens  is  seen 
to  be  heterogeneous,  and  in  many  cases  the  fibers  appear  to  be  in  contact. 
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DISCUSSION 


The  measured  critical  energy  release  rates  are  seen  to  have  a  definite 


crack  growth  rate  dependence.  Schapery  [12]  has  predicted  that  crack 
growth  rates  in  viscoelastic  material  may  be  described  by  an  equation  of 


the  form 


dL  = 
dt 


(1  +P 

C1GC 


(14) 


if  the  intrinsic  fracture  energy  and  strength  are  independent  of  crack 
speed.  The  exponent  depends  on  n,  which  is  a  material  constant  that 
describes  the  resin  creep  complaince  D(t)  through  an  equation  of  the 
form  D  =  Dg  +  D-j  tn.  The  value  of  n  =  0.1  is  obtained  from  Figure  2 
and  Equation  (14).  In  order  to  check  this  value,  creep/recovery  tests 
were  run  on  90°  (fiber  angle)  specimens  of  Scotchply  to  determine  the 
time  dependence  of  the  matrix.  An  exponent  of  n  =  0.05  -  0.10  was 
determined  for  these  tests,  which  is  consistent  with  the  fracture 
results  for  delamination  fracture  toughness. 

Delamination  fracture  toughness  of  Scotchply  has  been  measured  by 
Wu  [13]  using  a  center-notched  specimen  with  90°  fibers.  He  obtained  a 
value  of  1.75  ksi/Tn  for  Kc-  These  results  may  be  compared  to  our  results 
if  our  measure  dc  values  are  expressed  in  terms  of  Kc  using  the  relation¬ 
ship  1 

-  7 

K,c  ■  <EG> 

where  E  is  the  effective  modulus  and  has  been  shown  to  be  approximately 
equal  to  one  half  the  transverse  modulus  E22  for  an  orthotropic  material  [12]. 


Utilizing  the  manufacturer's  predicted  value  for  E22  of  1.4  x  10^  psi 
and  the  observed  range  of  Gc  values  of  3.0  -  5.5  lbs/in.  gives  a 
predicted  range  of  values  for  Kc  of  1.5  -  2.0  ksi/Tn,  in  good  agreement 
with  Wu's  results.  These  values  may  also  be  compared  with  values  measured 
for  metals  and  other  plastics,  as  summarized  in  Table  4. 

Several  assumptions  were  made  implicitly  in  the  analysis  that  appear 
to  be  justified  by  the  good  agreement  between  the  predictions  and 
measurements  for  the  nonlinear  beam  analysis,  as  seen  in  Figure  9.  These 
assumptions  include  plane  stress  and  "rigid  wall"  behavior  at  beam  end 
point,  which  in  our  specimen  is  the  crack  front.  The  actual  behavior 
in  this  region  may  only  be  described  by  a  more  complex  three  dimensional 
analysis  using  finite  elements;  however,  the  departure  from  rigid  wall 
behavior  appears  to  have  a  sufficiently  small  effect  on  the  overall 
strain  energy  stored  in  that  it  may  effectively  be  neglected.  Fiber 
breakage  in  the  other  fibers  could  also  have  introduced  an  error  into  the 
strain  energy  calculations  as  well  as  the  measured  critical  energy  release 
rate.  Simple  calculations  indicate  that  the  axial  stress  in  the  beam 
never  reached  the  failure  value;  furthermore,  almost  no  fiber  fracture 
was  noted  in  the  tests.  Finally,  as  a  result  of  transverse  strains 
(normal  to  the  loading  direction),  the  beam  has  a  compound  curvature 
which  would  tend  to  stiffen  the  beam  and  depend  on  beam  length.  Again, 
in  view  of  the  agreement  in  Figure  9,  this  effect  was  apparently  neglible 
for  the  specimens  used. 
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CONCLUSION 

A  new  approach  to  the  determination  of  delamination  fracture 
toughness  in  the  opening  mode  has  been  developed  using  a  double 
cantilever  beam  with  a  nonlinear  beam  analysis.  The  analysis  has  been 
confirmed  with  experimental  measurements  on  Scotchply  using  specimens 
of  three  thicknesses,  tested  over  a  wide  range  of  crack  growth  rates. 

The  measured  critical  energy  release  rate  of  3.0  -  5.5  Ibs/in  for  the 
range  of  crack  growth  rates  studies  has  been  analyzed  using  viscoelastic 
crack  growth  rate  theory  and  found  to  give  results  consistent  with  the 
viscoelastic  behavior  of  the  resin  as  determined  from  creep/recovery 
tests. 
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Table  1.  Nonl inear  Variables 
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3.8606 

14.9040 

.8472 

2.2416 
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oo 

00 
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00 
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Table  2.  Experimentally  Measured  Values  of  Crack  Length  (L),  Load  (P) 
and  Grip  Displacement  (?a)  for  Twelve  Ply  Sepcimen.  Also 
Tabulated  are  Various  Intermediate  Terms  Calculated  in 
Determining  G  from  the  Experimentally  Measured  Terms. 


L 

(in) 

2a 

(in) 

(lbs) 

A/L 

PL  /El 

WL 

El 

S 

Gc 

lbs/ in 

3.5 

2.22 

6.07 

.317 

1.0681 

.1596 

2.1238 

5.65 

4.0 

2.85 

5.48 

.3563 

1.2307 

.2045 

2.1560 

5.73 

4.5 

3.61 

4.97 

.4011 

1.4446 

.2644 

2.1919 

5.79 

5.0 

4.37 

4.62 

.4370 

1.6420 

.3198 

2.2338 

5.80 

5.5 

5.14 

4.25 

.4672 

1.8270 

.3719 

2.2824 

5.66 

6.0 

5.96 

3.99 

.4967 

2.0259 

.4285 

2.3410 

5.62 

6.5 

6.92 

3.81 

.5323 

2.2965 

.5051 

2.4231 

5.72 

7.0 

7.80 

3.67 

.5571 

2.5116 

.5647 

2.4866 

5.73 

7.5 

8.83 

3.54 

.5887 

2.8296 

.6491 

2.5761 

5.79 

8.0 

9.71 

3.33 

.6069 

3.0414 

.7027 

2.6342 

5.58 

8.5 

10.58 

3.17 

.6224 

3.2418 

.7515 

2.6892 

5.41 

9.0 

11.38 

2.91 

.6322 

3.3794 

.7840 

2.7272 

5.02 

9.5 

12.64 

2.97 

.6653 

3.9188 

.9044 

2.8809 

5.30 

10.0 

13.81 

3.01 

.6905 

4.4269 

1.0089 

3.0321 

5.52 
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Table  3.  Comparison  of  Measured  Values  of  Crack  Length  ( L)  To  Value 
Calculated  Using  Nonlinear  Beam  Analysis  and  Measured  Values 
of  Grip  Displacement  (2a)  and  Load  (P). 


Measured  Crack  Length  (in.)  Calculated  Crack  Length  (in.) 


3.50 

3.51 

4.00 

3.99 

4.50 

4.52 

5.00 

5.00 

5.50 

5.50 

6.00 

5.99 

6.50 

6.51 

7.00 

6.97 

7.50 

7.50 

8.00 

8.01 

8.50 

8.49 

9.00 

9.02 

9.50 

9.55 

10.00 


10.06 
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Table  4.  Typical  Fracture  Toughness  Parameters  (K^  and  Gc) 


Material 

KIC 

(ksi/Tn) 

Gc 

(lbs/in) 

Metals  [5,  101 

T i -6A1 -4V 

105 

689 

7075-T651 

22 

47 

4340 

55  -  90 

104  -  279 

2024  T3 

40 

152 

Thermoplastic  Materials  [ll,  12] 

Polymethyl  Methacrylate 

1.46  -  1.74 

6.05 

Polystyrene 

.894  -  1.00 

2.0 

Polyvinyl  Chloride 

1.45  -  2.05 

7.024 

Nylon  -  6,  6 

.462  -  .756 

1.428 

Polyethylene 

.756  -  1.069 

28.5 

Figure  la.  Wang/Mandell  specimen  for  delamination 
fracture  toughness  studies 


Figure  lb.  Double  cantilever  beam  specimen  used  in 
this  work  for  del  ami  nation  fracture 
toughness  studies 


A/L 


2 

Figure  3.  (PL  /El)  vs.  (a/L)  as  determined  from  nonlinear  beam  analysis 
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Figure  8.  Double  cantilever  beam  test  sequence. 


(a)  8  ply  specimen,  lOOx. 


(c)  12  ply  specimen,  lOOx. 


(d)  16  ply  specimens,  lOOx. 


Flgureslla-lld  (continued). 


